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"a^^ 

O^ i We review the formalism of the effective average action in quantum field theory 

O^ ' which corresponds to a coarse grained free energy in statistical mechanics. The 

associated exact renormalization group equation and possible nonperturbative ap- 
proximations for its solution are discussed. We describe in detail 0(Af)— symmetric 
scalar theories in two and three dimensions. These ideas are also applied to QCD 
where one observes the consecutive emergence of mesonic bound states and sponta- 
neous chiral symmetry breaking as the coarse graining scale is lowered. We finally 
present a study of the chiral phase transition in two flavor QCD. A precision es- 
timate of the universal critical equation of state for the three— dimensional 0(4) 
Heisenberg model is presented. We explicitly connect the 0(4) universal behavior 
near the critical temperature and zero quark mass with the physics at zero tem- 
perature and a realistic pion mass. For realistic quark masses the pion correlation 
length near T^ turns out to be smaller than its zero temperature value. 
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1 Effective average action 
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Qs , Quantum chromodynamics (QCD) describes qualitatively different physics at 

Q> ■ different length scales. At short distances the relevant degrees of freedom are 

quarks and gluons which can be treated perturbatively. At long distances we 
O ' observe hadrons, and an essential part of the dynamics can be encoded in 

the masses and interactions of mesons. Any attempt to deal with this situa- 
tion analytically and to predict the meson properties from the short distance 
(-* , physics (as functions of the strong gauge coupling as and the current quark 

masses ruq) has to bridge the gap between two qualitatively different effective 
descriptions. Two basic problems have to be mastered for an extrapolation 
from short distance QCD to mesonic length scales: 

h '. 

Cu ^ • The effective couplings change with scale. This does not only concern the 

running gauge coupling, but also the coefficients of non~renormalizable 
operators as, for example, four quark operators. Typically, these non- 
renormalizable terms become important in the momentum range where 
a^ is strong and deviate substantially from their perturbative values. 
Consider the four-point function which results after integrating out the 
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gluons. For heavy quarks it contains the information about the shape 
of the heavy quark potential whereas for hght quarks the comphcated 
spectrum of hght mesons and chiral symmetry breaking are encoded in 
it. At distance scales around Ifm one expects that the effective action 
resembles very little the form of the classical QCD action which is relevant 
at short distances. 

• Not only the couplings, but even the relevant variables or degrees of 
freedom are different for long distance and short distance QCD. It seems 
forbiddingly difficult to describe the low-energy scattering of two mesons 
in a language of quarks and gluons only. An appropriate analytical field 
theoretical method should be capable of introducing field variables for 
composite objects such as mesons. 

A conceptually very appealing idea for our task is the block-spin actionB'cl It 
realizes that physics with a given characteristic length scale I is conveniently 
described by a functional integral with an ultraviolet (UV) cutoff A for the 
momenta. Here A should be larger than l^^ but not necessarily by a large 
factor. The Wilsonian effective action 5^ replaces then the classical action 
in the functional integral. It is obtained by integrating out the flucf. 
with momenta q^ ^ A^. An exact renormalization group equatior 
describes how S^ changes with the UV cutoff A. 

We will usahere the somewhat different but related concept of the effective 
average actionQ Tk which, in the language of statistical physics, is a coarse 
grained free energy with coarse graining scale k. The effective average action 
is based on the quantum field theoretical concept of the effective actionEl F 
which is obtained by integrating out all quantum fluctuations. The effective 
action contains all information about masses, couplings, form factors and so 
on, since it is the generating functional of the IPI Green functions. The 
field equations derived from F are exact including all quantum effects. For 
a field theoretical description of thermal equilibrium this concept is easily 
generalized to a temperature dependent effective action which includes now 
also the thermal fluctuations. In statistical physics F describes the free energy 
as a functional of some convenient (space dependent) order parameter, for 
instance the magnetization. In particular, the behavior of F for a constant 
order parameter (the effective potential) specifies the equation of state. The 
effective average action Ffc is a simple generalization of the effective action, 
with the distinction that only quantum fluctuations with momenta q^ <; k^ 
are included. This can be achieved by introducing an explicit infrared cutoff 
~ fc in the functional integral defining the partition function (or the generating 
functional for the n-point functions). Typically, this IR-cutoff is quadratic in 
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the fields and modifies the inverse propagator, for example by adding a mass- 
like term ~ k^. The effective average action can then be defined in complete 
analogy to the effective action (via a Legendre transformation of the logarithm 
of the partition function) . The mass-like term in the propagator suppresses 
the contributions from the small momentum modes with q^ ^ k'^ and Ffc 
accounts effectively only for the fluctuations with q^ ^ k"^. 

Following the behavior of Ffe for different k is like looking at the world 
through a microscope with variable resolution: For large k one has a very 
precise resolution ~ k~^ but one also studies effectively only a small volume 
^ k"^. Taking in QCD the coarse graining scale k much larger than the 
confinement scale guarantees that the complicated nonperturbative physics 
does not play a role yet. In this case, F^ will look similar to the classical 
action, typically with a running gauge coupling evaluated at the scale k. (This 
does not hold for Green functions with much larger momenta p^ 3> k^ where 
the relevant IR cutoff is p, and the effective coupling is a,s{p).) For lower k the 
resolution is smeared out and the detailed information of the short distance 
physics can be lost. (Again, this does not concern Green functions at high 
momenta.) On the other hand, the "observable volume" is increased and long 
distance aspects such as collective phenomena become visible. In a theory 
with a physical UV cutoff A we may associate Fa with the classical action 
S since no fluctuations are effectively included. By deflnition, the effective 
average action equals the effective action for A: = 0, Fq = F, since the infrared 
cutoff is absent. Thus F^ interpolates between the classical action S and the 
effective action F as /c is lowered from A to zero. The ability to follow the 
evolution A: ^ is equivalent to the ability to solve the quantum field theory. 

For a formal description we will consider in the first two sections a model 
with real scalar fields x" (the index a labeling internal degrees of freedom) in 
d Euclidean dimensions with classical action ^[x]- We define the generating 
functional for the connected Green functions by 

Wk[J]=\nJvxcxpi-Sk[x]+ ld''xJa{x)x''{x)\ (1) 

where we have added to the classical action an IR cutoff A^S 

Sk[x]^S[x] + ASk[x] (2) 

which is quadratic in the fields and conveniently formulated in momentum 
space 



^SM^lj ^,Rk{q^)Xa{-q)x^ 



(q) ■ (3) 
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Here Ja are the usual scalar sources introduced to define generating functional 
and Rk{q^) denotes an appropriately chosen (see below) IR cutoff function. 
We require that Rk{q^) vanishes rapidly for q^ ^ k^ whereas for q^ ^ k^ it 
behaves as Rk{q^) — k^ ■ This means that all Fourier components of x" with 
momenta smaller than the IR cutoff k acquire an effective mass rrics ~ k and 
therefore decouple while the high momentum components of x°' should not 
be affected by Rk- The classical fields 

*"-«-^ (4) 

now depend on fc. In terms of Wk the effective average action is defined via a 
Legendre transformation 

Tfc [$] ^-Wk[J]+ f d'^xJa (x)^'' (x) - ASk m . (5) 

In order to define a reasonable coarse grained free energy we have subtracted 
in (P) the infrared cutoff piece. This guarantees that the only difference 
between Ffc and F is the effective IR cutoff in the fluctuations. Furthermore, 
this has the consequence that F^ does not need to be convex whereas a pure 
Legendre transform is always convex by definition. (The coarse grained free 
energy becomes convexB only for /c — j 0.) This is very important for the 
description of phase transitions, in particular first order ones. One notes 

lim Rk{q^) = => lim Ffc[$] = F[$] 

(6) 

lim Rk{q^) =oo^ lim Ffc[$] = S[^] 

fc— >A fc^A 

for a convenient choice of Rk like 

Rk{q)^ Z^M ——rj^r——rj]^r- ■ (7) 

Here Z^ k denotes the wave function renormalization to be specified below 
and we will often use for Rk the limit A — j oo 

i?.(.^) = J^- (8) 

We note that the property Fa = 5 is not essential since the short distance 
laws may be parameterized by Fa as well as by S. In addition, for momentum 
scales much smaller than A universality implies that the precise form of Fa is 
irrelevant, up to the values of a few relevant renormalized couplings. 

A few properties of the effective average action are worth mentioning: 
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1. All symmetries of the model which are respected by the IR cutoff A^S are 
automatically symmetries of Tk- In particular, this concerns translation 
and rotation invariance and one is not plagued by many of the problems 
encountered by a formulation of the block-spin action on a lattice. 

2. In consequence, Tk can be expanded in terms of invariants with respect 
to these symmetries with couplings depending on k. For the example of 
a scalar theory one may use a derivative expansion {p = <I>"<I>a/2) 

Tfc = Jd^x lukip) + lz^Ap)d''^ad^^" + ■■■] (9) 

and expand further in powers of p 

1- 2 1_ 3 

Uk{p) = ^>^k (p - Po{k)) + g7fc(p-Po(fc)) +■■■ 

z<i>,k{p) = z^Apo) + z'^Ap^) (p- Po) + ■■■ (10) 

where po denotes the (fc-dependent) minimum of the effective average 
potential Uk{p)- We see that Tk describes infinitely many running cou- 
plings. (Z$,fc in (Q) can be identified with Z,i,_k{po)-) 

3. There is no problem incorporating^chiral fcrmions since a chirally invari- 
ant cutoff Rk can be formulated.E3llil 

4. Gauge |— ■ j— Jlieaides, j— . |-eaii| ■— ■ be formulated along similar 



line£3'liia't3lla'li3ll3'E3tal£J'Eil even though A^S' may not be gauge 
invariant Jj In this case the usual Ward iidentities receive corrections for 
which one can derive closed expressions. t^ These corrections vanish for 
fc^O. 

5. The high momentum modes are very effectively integrated out because 
of the exponential decay of Rk for q^ ^ fc^. Nevertheless, it is some- 
times technically easier to use a cutoff without this fast decay property 
(e.g. Rk ^ k^ or Rk ~ k'^/q^). In the latter case one has to be careful 
with possible remnants of an incomplete integration of the short distance 
modes. Also our cutoff does not introduce any non-analytical behavior 
as would be the case for a sharp cutoff.L^ 

6. Despite a similar spirit and many analogies there remains also a con- 
ceptual difference to the Wilsonian effective action 5^. The Wilsonian 

"For a manifestly gauge invariant forinulatioB-Jn terms of Wilson loops seeB3. 
''For results using a sharp cutoff function sccr'1. 
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effective action describes a set of different actions (parameterized by A) 
for one and the same theory — the n-point functions are independent of 
A and have to be computed from S"^ by further functional integration. 
In contrast, Tk describes the effective action for different theories — for 
any value of k the effective average action is related to the generating 
functional of a theory with a different action Sk ~ S + AkS. The n— 
point functions depend on k. Th£_.Wilsonian effective action does not 
generate the IP/ Green functions.E3 

7. Because of the incorporation of an infrared, cutoff , Ffc is closely related 
to an effective action for averages of fields ,0 where the average is taken 
over a volume ^ k'^ . 



2 Exact renormalization group equation 

The dependence of Tk on the coarse graining scale k is governed by an exact 
renormalization group equation (ERGE)E£l 

dtV^m = i Tr I [ri^) [$] + i?fe] "' dtRK^ . (11) 

Here t — ln(/c/A) with some arbitrary momentum scale A, and the trace 
includes a momentum integration as well as a summation over internal indices, 
Tr = J , ?j Tlia- ^^^ second functional derivative F^ ' denotes the exact 
inverse propagator 

The flow equation (|ll|) can be derived from (|^) in a straightforward way using 

dtTkU = -dtWk\j-dAkSm 

= iTr{a,i?..P^f } (13) 



and 






5J"-{-q)5J''{q') 
5^Wk S^ (Tk + AkS) 



SJa{-q)SMq') 6M-Q')^'^ciq") 



SacSqq" ■ (14) 
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It has the form of a renormaUzation group improved one-loop expression^ 
Indeed, the one-loop formula for r^ reads 

Tfc [$] = 5[$] + i Tr In (s^^) [$] + i?,.) (15) 

with S'"^ the second functional derivative of the classical action, similar 
to (p^. (Remember that S''-^-' is the field dependent classical inverse propa- 
gator. Its first and second derivative with respect to the fields describe the 
classical three- and four-point vertices, respectively.) Taking a i-derivative 
of (ll^) gives a one-loop flow equation very similar to (O) with F^ ' replaced 



fc 
by S^'^K It may seem surprising, but it is nevertheless true, that the renormal- 

ization group improvement S^"^' -^ FJ. promotes the one-loop flow equation 
to an exact nonperturbative flow equation which includes the effects from all 
loops as well as all contributions which are non-analytical in the couplings like 
instantons, etc.! For practical computations it is actually often quite conve- 
nient to write the flow equation (O) as a formal derivative of a renormalization 
group improved one-loop expression 

dtrk^^TYdt\n(r[^'^+Rk) (16) 

with dt acting only on i?^ and not on F^ , i.e. dt = {dRk/dt) (d/dRk). Flow 
equations for n-point functions follow from appropriate functional derivatives 
of ( [ll|) or (|l^) with respect to the flelds. For their derivation it is sufficient to 
evaluate the corresponding one-loop expressions (with the vertices and prop- 
agators derived from F^) and then to take a formal 9* -derivative. (If the 
one-loop expression is finite or properly regularized the 9t -derivative can be 
taken after the evaluation of the trace.) This permits the use of (one-loop) 
Feynman diagrams and standard perturbative techniques in many circum- 
stances. Most importantly, it establishes a very direct connection between 
the solution of flow-equations and perturbation theory. If one uses on the 
right hand side of (O) a truncation for which the propagator and vertices ap- 



(2) 

pearing in F^ are replaced by the ones derived from the classical action, but 
with running fc-dependent couplings, and then expands the result to lowest 
non-trivial order in the coupling constants one recovers standard renormaliza- 
tion group improved one-loop perturbation theory. The formal solution of the 
flow equation can also be employed for the development of a systematically 
resummed perturbation theory.EZI 

For a choice of the cutoff function similar to (pi) the momentum integral 
contained in the trace on the right hand side of the flow equation is both 
infrared and ultraviolet finite. Infrared finiteness arises through the presence 
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of the infrared regulator ^ Rk- We note that all eigenvalues of the matrix 
T]^' + Rk must be positive semi-definite. The proof follows from the observa- 
tion that the functional Tk + A^S* is convex since it is obtained from Wk by 
a Legendre transform. On the other hand, ultraviolet finiteness is related to 
the fast decay oi dtRk for q^ ^ k^. This expresses the fact that only a narrow 
range of fluctuations with q^ ~ k^ contributes effectively if the infrared cutoff 
k is lowered by a small amount. If for some other choice of Rk the right hand 
side of the flow equation would not remain UV finite this would indicate that 
the high momentum modes have not yet been integrated out completely in 
the computation of Tk- Since the fiow equation is manifestly finite this can 
be used to define a regularization scheme. The "ERGE-scheme" is specified 
by the fiow equation, the choice of Rk and the "initial condition" Fa. This is 
particularly important for gauge theories where other regularizations in four 
dimensions and in the presence of chiral fermions are difficult to construct. 
For gauge theories Fa has to obey appropriately modified Ward identities. In 
the context of perturbation theory a first proposal. hDw to regularize gauge 
theories by use of flow equations can be found in.EJ'tj We note that in con- 
trast to previous versions of exact renormalization group equations there is no 
need in the present formulation to construct an ultraviolet momentum cutoff 
— a task known to be very difficult in non-Abelian gauge theories. 

Despite the conceptual differences between the Wilsonian effective action 
S'a' and the effective average action F^ the exact flow equations describing 
the A-dependence of SY and the fc-dependcnce of Fj, are simply related. 
Polchinski's continuum version of the Wilsonian flow equationd can be trans- 
formed into^Ul]) by means of a Legendre transform and a suitable variable 
redefinition.^ 

Even though intuitively simple, the replacement of the (RG-improved) 
classical propagator by the full propagator turns the solution of the fiow equa- 
tion ( [ll|) into a difficult mathematical problem: The evolution equation is a 
functional differential equation. Once Ffc is expanded in terms of invariants 
(e.g. Eqs.(||), (p^) this is equivalent to a coupled system of non-hnear partial 
differential equations for infinitely many couplings. General methods for the 
solution of functional differential equations are not developed very far. They 
are mainly restricted to iterative procedures that can be applied once some 
small expansion parameter is identified. This covers usual perturbation the- 
ory in the case of a small coupling, the 1/A'^-expansion or expansions in the 
dimensionality A — d or 2 — d. It may also be extended to less familiar expan- 
sions like a derivative expansion which is related in critical three dimensional 
scalar theories to a small anomalous dimension. In the absence of a clearly 
identified small parameter one nevertheless needs to truncate the most gen- 
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eral form of r^ in order to reduce the infinite system of coupled differential 
equations to a (numerically) manageable size. This truncation is crucial. It is 
at this level that approximations have to be made and, as for all nonperturba- 
tive analytical methods, they are often not easy to control. The challenge for 
nonperturbative systems like low momentum QCD is to find flow equations 
which (a) incorporate all the relevant dynamics such that neglected effects 
make only small changes, and (b) remain of manageable size. The difficulty 
with the first task is a reliable estimate of the error. For the second task the 
main limitation is a practical restriction for numerical solutions of differential 
equations to functions depending only on a small number of variables. The 
existence of an exact functional differential flow equation is a very useful start- 
ing point and guide for this task. At this point the precise form of the exact 
flow equation is quite important. Furthermore, it can be used for systematic 
expansions through enlargement of the truncation and for an error estimate in 
this way. Nevertheless, this is not all. Usually, physical insight into a model 
is necessary to device a useful nonperturbative truncation! 

So far, two complementary approaches to nonperturbative truncations 
have been explored: an expansion of the effective Lagrangian in powers of 
derivatives (p = i<I)a$") 

(17) 
or one in powers of the fields 

[$(x,)-*o] 1 r^"^(xi,...,x„). (18) 




If one choosesl^ $0 as the fc-dependent VEV of $, the series ( |18|) starts 



in) 

effectively at ri = 2. The flow equations for the IP I n-point functions F), 



are obtained by functional differentiation of (|1 1|) . Such flow- equations have 
been discussed earlier from a somewhat different viewpoint D They_can also 
be interpreted as a differential form of Schwinger-Dyson equations .l2I 

The formation of mesonic bound states, which typically appear as poles in 
the (Minkowskian) four-quark Green function, is most efficiently described by 
expansions like (|l8|). This is also the form needed to compute the nonpertur- 
bative mcttijipit|mn dependence of the gluon propagator and the heavy quark 
potential.cJ'llII'ca On the other hand, a parameterization of F^ as in (|lj) seems 
particularly suited for the study of phase transitions. The evolution equation 



'^See also B3 for the importance of expanding around # = <J>o instead of <!> = 0. 
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for the average potential Uk follows by evaluating ( [l7| ) for constant 4>. In t 
limit where the ^-dependence of Z^k is neglected and Y$,fc ~ one find 
for the 0(iV)-synimetric scalar model 

„ ^^ , , I f dU dRk( N-\ 



2 J (27r)'i dt \Z^,kq^ + Rk + U'^ 



k 



^ Z^M^ + Ru {q) + UI + 2pU'l ) ^^^^ 

with C/^ = ^^, etc. One observes the appearance of p-dependent mass 
terms in the effective proaagators of the right hand side of ( p^ ) . Once 77$ = 
—dt In Z^^k is determinecfl in terms of the couplings parameterizing Uk this 
is a partial differential equation for a funptianiiwdepending on two variables 
k and p which can be solved numerically.OOm'LJ (The Wilson-Fisher fixed 
point relevant fcr-a, second order phase transition [d — 3) corresponds to a 
scaling solutionEj' L3 where dtUk ~ 0) A suitable truncation of a flow equation 
of the type ( |l7|) will play a central role in the description of chiral symmetry 
breaking below. 

It should be mentioned at this point that the weakest point in the ERGE 
approach seems to be a reliable estimate of the truncation error in a non- 
perturbative context. This problem is common to all known analytical ap- 
proaches to nonperturbative phenomena and appears often even within sys- 
tematic (perturbative) expansions. One may hope that the existence of an 
exact flow equation could also be of some help for error estimates. An obvious 
possibility to test a given truncation is its enlargement to include more vari- 
ables — for example, going one step higher in a derivative expansion. This 
is similar to computing higher orders in perturbation theory and limited by 
practical considerations. As an alternative, one may employ different trun- 
cations of comparable size — for instance, by using different definitions of 
retained couplings. A comparison of the results can give a reasonable picture 
of the uncertainty if the used set of truncations is wide enough. In this con- 
text we should also note the dependence of the results on the choice of the 
cutoff function Rk{q)- Of course, for fc -^ the physics should not depend 
on a particular choice of Rk and, in fact, it does not for full solutions of (pj]). 
Different choices of Rk just correspond to different trajectories in the space 
of effective average actions along which the unique IR limit r[$] is reached. 
Once approximations are used to solve the ERGE (pi), however, not only the 
trajectory but also its end point will depend on the precise definition of the 
function Rk. This is very similar to the renormalization scheme dependence 
usually encountered in perturbative computations of Green functions. One 
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may use this scheme dependence as a tool to study the robustness of a given 
approximation scheme. 

Before applying a new nonperturbative method to a complicated theory 
like QCD it should be tested for simpler models. A good criterion for the ca- 
pability of the ERGE to deal with nonperturbative phenomena concerns the 
critical behavior in three dimensional scalar theories. In a first step the well 
known results of other methods for the critical exponents have been repro- 
duced within a few percent accuracy.t3 The ability of the method to produce 
new results has been demonstrated by thejiomputation of the critical equation 
of state for Ising and Heisenberg modelS which has been verified by lattice 
simulatiDns.LO This has been extended to first order transitions in matrix 
modelsa or for the Abelian Higgs model relevant for superconductors.HE^ 
Analytical investigations of the high temperature phase transitions in d = 4 
scalar theories (0(7V)pmodels) have correctly described the second order na- 
ture of the transition,LJ in contrast to earlier attempts within high tempera- 
ture perturbation theory. 

For an extension of the flnw pmia.tinns to A helian and non- Abelian gauge 
theories we refer the reader to.E3'L3EJ't§llj'OEfl The other necessary piece 
for a description of low-energy QCD, namely the transition from fundamen- 
tal (quark and gluon) degrees of freedom to composite (meson) fields within 
the framework of the ERGE can be found in.E3 We will describe the most 
important aspects of this formalism for mesons below. 

3 Scalar 0{N) models in two and three dimensions 

The universal critical behavior of many systems of statistical mechanics is 
described by the field theory for scalars with 0{N) symmetry. This covers 
the gas-liquid and many chemical transitions described by Ising models with 
a discrete symmetry Z2 = 0(1), superfluids with continuous abelian symme- 
try 0(2), Heisenberg models for magnets with A^ = 3, etc. In 2 < d < 4 
dimensions all these models have a continuous second order phase transition. 
In two dimensions one observes a second, order transition for the Ising model, 
a Kosterlitz-Thouless phase transitioned for N ~J2. and no phase transition 
for non-abelian symmetries A^ > 3. It is knownH that a continuous sym- 
metry cannot be broken in two dimensions in the sense that the expectation 
value of the unrenormalized scalar field vanishes in the limit of vanishing 
sources, $a = (xaix)) — 0. We will see that all this diversity can already 
be described using an extremely simple truncation of the flow equation for 
0(A'^)-symmetric linear $^-model. 



We start from the flow equation (19) for arbitrary d and introduce the 
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anomalous dimension 



and the threshold function 



d 
m = -^lnZ$,fc (20) 



itiw; ,,) = '-v-/k-^z-^i I .:i, ,, _rn. . ,. (21) 



where 



^-1 ^ 2'^+\'^/^r{d/2) . (22) 



The threshold function decays fast for w ^ 1 and therefore expresses the 
decoupling of heavy particles with mass Af j 3> Z^^kk^ from the renormaliza- 
tion group flow. Its precise shape depends on the choice of Rk. It contains 
a term Hnear in r/$ from Z^^f,dtRk = dt{Rk/Z^,k) - m{Rk/Z^^k), cf. eq. (§). 
In terms of this threshold function one has 



Z^ kk Z<^ kk 



dtUk = 2vak''\ it{ C.,'^^,' ;m) + iN- ir^i^i^;^^) . (23) 



For given 77$ (fc) (see below) this is a partial differential equation for a function 
of two variables Uk{p)- It may be solved numerically for arbitrary d and N . 
For given "initial conditions" specifying Ua{p) for some microscopic length 
scale A~^ the free energy results from the solution for A; ^ 0. In particular, 
the expectation value (%) corresponds to the minimum of Uk for fc — > 0. De- 
noting the running location of the potential minimum by po{k), spontaneous 
symmetry breaking occurs for po(0) > 0, whereas po{0) = corresponds to the 
symmetric or disordered phase. The dependence on temperature T is reflected 
by the T-dependence of the mass term or the location of the minimum at the 
scale A. Varying po{A) one may detect the different phases. In particular, 
a second order phase transition corresponds to a critical value po{A) = po* 
such that for po(A) > po* spontaneous symmetry breaking occurs, whereas for 
po(A) < po* one ends in the symmetric phase with po(0) = and Uq{0) > 0. 
Near the critical temperature Tc one may linearize po(A) = pa* + A{Tc — T). 
Instead of describing the results of a numerical solutionl23 of eq. ( p3| ) we 
make here a very simple polynomial approximation to [7^, namely 

Uk{p)^\\k{p-po{k)f . (24) 

All chaxacteristic features of the flow can already be seen in the approxi- 
mationfl Furthermore, we also approximate ^0(10; 77$) by l^iw) = Iq{'W]0) 
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which is well justified for small ?/$. The flow of the potential niinimum can 
be inferred from the identity 

= jUlip^ik)) = dtUl{p^{k)) + C/^'(po(fc))9tPo(fc) . (25) 

Here dtUl{p) is the partial i-derivative with p held fixed which is computed 
by differentiating eq. ( p3| ) with respect to p. Defining the additional threshold 
functions by 

d 

Zf,+i(u;;77$) = -i— /^(u;;r,$), n>l (26) 

n ow 

and constants Zq = /jj(0; 0) one obtains 

d,p, = 2v,k'-^Z^^,Llti^^) + (iV - 1);A . (27) 

We conclude that po{k) always decreases as the infrared cutoff k is lowered. 
For d = 2 and N > 2 only a Z-factor increasing without bounds can prevent 
po{k) from reaching zero at some value A; > 0. We will see that this happens 
for TV = 2 in the low temperature phase. 

It is convenient to introduce renormalized dimensionless fields and a di- 
mensionless potential 

p^Zg.,kk^^''p, u^k^^Uk (28) 

and corresponding renormalized dimensionless couplings as 

K = Z$.fefc2-rfpg ^ X^Z^lk'^-'^X. (29) 

This yields the scaling formE3 of the flow equation (23) 
dtu\- = —du + (d — 2 + r/^)pu' 

+ 2vJ li{u' + 2pu"- m) + iN- l)l^iu', rj^) I (30) 

where u' denotes du/ dp. All explicit scale and wave function renormalization 
factors have disappeared from this partial differential equation^ for u{p,t). 
For our simple truncation one has 

w=1a(p-k)2 (31) 



"^The critical soluti«n— (9tu = is an ordinary differential equation for u{p) which can be 



solved numericalljo'l 
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and the flow equations for k and A read 

dtK = P^ = {2-d~ r]^)K + 2vd< 3lf{2XK) + {N - l)lf \ (32) 

5tA = /?A = (d - 4 + 277<i.)A + 2vdX^l 9^2 (2Ak) + {N - 1)4 \ ■ (33) 

In this truncation the anomalous dimension ?/$ is given byij 

r;<(. = l^A2Km^_2(0,2AK) (34) 

where TO2 2 is another threshold function with the property 

7?$ = for d = 2 , Ak > 1 . (35) 

Attk 

We want to demonstrate next that the two differential equations (|3^), 
( |33| ) describe all phase transitions in two or three dimensions correctly. For 
d = 3 one finds a fixed point (k^,A*) where [3^. = P\ = 0. From the generic 
form fix = —A + A^(ci + C2(Ak)) one concludes that A essentially corresponds 
to an infrared stable coupling which is attracted towards its fixed point value 
A* as k is lowered. On the other hand, /J^ — —n + C3 + C4(Ak) shows that 
K is essentially an infrared unstable or relevant coupling. Starting for given 
Aa with KA = '«*(Aa) + Ska, Ska = a{Tc — T), a > 0, one either ends in 
the symmetric phase for Ska < or spontaneous symmetry breaking oc- 
curs for Ska > 0. The fixed point or, more generally, the scaling solution 
with dtu — corresponds precisely to the critical temperature of a second 
order phase transition. Critical exponents can be computed from solutions 
in the vicinity of the scaling solution. The index v characterizes the diver- 
gence of the correlation length for T -^ T^, i.e., ^ ^ m'^^ ^ \T ~ Tc\~'^ 
with m|j — Hmfe^o Z^]. [[/^(po) + 2pot^fc (Po)]- It corresponds to the negative 
eigenvalue of the "stability matrix" A^ = {dj3i/dXj){K^,,X^) with A^ = (k. A). 
(This can be generalized for more than two couplings.) The critical expo- 
nent 77 determines the long distance behavior of the two-point function for 
T = Tc- It is simply given by the anomalous dimension at the fixed point, 
77 = 77$ (k*. A*). It is remarkable that already in a very simple poljppmial 
trunmtion the critical exponents come out with reasonable accuracy£j (See 
alsocS for a discussion of the iV = 1 case in three dimensions.) 

In two dimensions the term linear in k vanishes in (3,^. This changes the 



faro: submitted to World Scientific oti December 18, 1998 14 



fixed point structure dramatically as can be seen from 

lim /?« = ^^ (36) 

where l\ — 1 was used. Since /3k is always positive for k = a fixed point 
requir£S| itbat [3^ becomes negative for large k. This is the case for the Ising 
modelE3'E3 where A^ = 1. On the other hand, for a non-abelian symmetry 
with A^ > 3 no fixed point and therefore no phase transition occurs. The 
location of the minimum always reaches zero for some value fcs > 0. The only 
phase corresponds to a linear realization of 0{N) with N degenerate masses 
mR ~ fcs- It is interesting to note that the limit k — > oo describes the non- 
linear sigma model. The non-abelian coupling g of the non-linear model is 
related to k by g^ — 1/(2k) and eq. ( p6| ) reproduces the standard one-loop 
beta function for g 



dg^ _ N^2 



/ (37) 



which is characterized by asymptotic freedom.ES The "confinement scale" 
where the coupling g becomes strong can be associated with kg. The strongly 
interacting physics of the non-linear model finds a simple description in terms 
of the symmetric phase of the linear 0(Af)-model!Ll 

Particularly interesting is the abelian continuous symmetry for N — 2. 
Here P^ vanishes for k — > oo and k becomes a marginal coupling. Usingj-the 
exact renormalization group approach outlined above one actually findalJ a 
behavior consistent with a second order phase transition with rj ~ 0.25 near 
the critical trajectory. The low temperature phase {k\ > k*) is special since 
it has many characteristics of the phase with spontaneous symmetry break- 
ing, despite the fact that poik — > 0) must vanish according to the Mermin- 
Wagner theorem.E^ There is a massless Goldstone-type boson (infinite corre- 
lation length) and one massive mode. Furthermore, the exponent rj depends 
on KA or the temperature (cf. eq. (|3g)), since k flows only marginally. JThese 
are the characteristic features of a Kosterlitz-Thouless phase transition.cil The 
puzzle of the Goldstone boson in the low temperature phase despite the ab- 
sence of spontaneous symmetry breaking is solved by the observation that the 
wave function renormalization never stops running with k: 

1/2 

Even though the renormalized field xb^ = ■^^ fcX acquires a non-zero expec- 
tation value (xr) — \/2k, for A: ^ the unrenormalized order parameter 
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vanishes due to the divergence of Z$_fc, 

«'»^yf(x)*- w 

Also the inverse Goldstone boson propagator behaves as (g2)i-i/(87r«;) g^j^^j 
circumvents Coleman's no-go theoremE3 for free massless scalar fields in two 
dimensions, ft is remarkable that all these features arise from the solution of a 
simple one-loop type equation without ever invoking nonperturbative vortex 
configurations. This enhances our confidence that the nonperturbative aspects 
are indeed already caught by relatively simple truncations of the effective 
average action. 

Having found all important qualitative features of the phase transition al- 
ready in a very rough truncation the way is open for systematic quantitative 
improvements. We report here only the results of a recent investigation re- 
garding the critical behavior in the universality class of the three-dimensional 
fsing model in next to leading order of a systematic derivative expansionB3 
This truncation is given by eq. (H^) with Y^,k = and used for a numeri- 
cal solution of the flow equation. One expects the flow of the potential and 
the exponent v to be described quite accurately in this approximation: Only 
the momentum dependence of the wave function renormalization is neglected 
which is governed by the small anomalous dimension. The uncertainties for 
T] are more substantial. One finds v — 0.631, 77 = 0.047 to be compared with 
the averages of several other methods v — 0.630, 77 = 0.035 which are pre- 
sumably more precise. Results of a similar quality are found for the universal 
amplitude ratios and the Widom scaling form of the critical equation of state. 

4 Chiral symmetry breaking in QCD 

The strong interaction dynamics of quarks and gluons at short distances or 
high energies is successfully described by quantum chropipdynamics (QCD). 
One of its most striking features is asymptotic freedomEj which makes per- 
turbative calculations reliable in the high energy regime. On the other hand, 
at scales around a few hundred MeV confinement sets in. As a consequence, 
the low-energy degrees of freedom in strong interaction physics are mesons, 
baryons and glueballs rather than quarks and gluons. When constructing 
effective models for these IR degrees of freedom one usually relies on the 
symmetries of QCD as a guiding principle, since a direct derivation of such 
models from QCD is still missing. The most important symmetry of QCD, 
its local color SU{3) invariance, is of not much help here, since the IR spec- 
trum appears to be color neutral. When dealing with bound states involv- 
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ing heavy quarks the so called "heavy quark symmetry" may be invoked to 
obtain approximate symmetry relations between IR observables.til We will 
rather focus here on the light scalar and pseudoscalar meson spectrum and 
therefore consider QCD with only the light quark flavors u, d and s. To 
a good approximation the masses of these three flavors can be considered 
as small in comparison with other typical strong interaction scales. One 
may therefore consider the chiral limit of QCD (vanishing current quark 
masses) in which the classical QCD Lagrangian does not couple left- and 
right-handed quarks. It therefore exhibits a global chiral invariance under 
ULiN) X UniN) = SUl{N) x SUr{N) x Uv{l) x Ua{1) where N denotes the 
number of massless quarks {N = 2 or 3) which transform as 

i^R=^-^^P^URi^n; Ur(.Ur{N) 

^I;^ = 1±J1^-^Ul^Pl\ Ul^Ul{N). (40) 

Even for vanishing quark masses only the diagonal SUv{N) vector-like sub- 
group can be observed in the hadron spectrum ("eightfold way"). The sym- 
metry SUl{N) X SUr{N) must therefore be spontaneously broken to SUv{N) 

SUl{N) X SUr{N) — ^ SUl+r{N) = SUv{N) . (41) 

Chiral symmetry breaking is one of the most prominent features of strong 
interaction dynamics and phenomenologically well established,E3 though a 
rigorous derivation of this phenomenon starting from first principles is still 
missing. In particular, the chiral symmetry breaking ( [4l| ) predicts for TV = 3 
the existence of eight light parity-odd (pseudo-)Goldstone bosons: 7r°, n^, 

K'^ , K , K and 77. Their comparably small masses are a consequence of 
the explicit chiral symmetry breaking due to small but non-vanishing cur- 
rent quark masses. The axial Abelian subgroup Ua{^) = Ul-r{1) is bro- 
ken in the quantum theory by an anomaly of the axial- vector iiurrent. This 
breaking proceeds without the occurrence of a Goldstone boson.E^ Finally, the 
Uv{^) — Ul+r{1) subgroup corresponds to baryon number conservation. 
The light pseudoscalar and scalar mesons are thought of as color neutral 

quark-antiquark bound states ^"^ ~ i'L'^R: a,b = 1, . . . ,N, which therefore 



transform under chiral rotations (|40D as 

$ — > Ur^uI . (42) 

Hence, the chiral symmetry breaking pattern (M^) is realized if the meson 
potential develops a VEV 

{<i>-'')=aoS-'>- ao^O. (43) 
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One of the most crucial and yet unsolved problems of strong interaction dy- 
namics is to derive an effective field theory for the mesonic degrees of freedom 
directly from QCD which exhibits this behavior. 

5 A semi quantitative picture 

Before turning to a quantitative description of chiral symmetry breaking using 
flow equations it is perhaps helpful to give a brief overview of the relevant 
scales which appear in relation to this phenomenon and the physical degrees 
of freedom associated to them. Some of this will be explained in more detail in 
the remainder of these lectures whereas other parts are rather well established 
features of strong interaction physics. 

At scales above approximately 1.5 GeV, the relevant degrees of freedom 
of strong interactions are quarks and gluons and their dynamics appears to be 
well described by perturbative QCD. At somewhat lower energies this changes 
dramatically. Quark and gluon bound states form and confinement sets in. 
Concentrating on the physics of scalar and pseudoscalar mesons |j there are 
three important momentum scales which appear to be rather well separated: 

• The compositeness scale fc$ at which mesonic tpip bound states form 
because of the increasing strength of the strong interaction. It will turn 
out to be somewhere in the range (600 — 700) MeV. 



• The chiral symmetry breaking scale k^sB at which the chiral conden- 
sate (ipTp°'\ or <^<i)°^^ assumes a non-vanishing value, therefore breaking 

chiral symmetry according to (|4l|). This scale is found to be around 
(400 — 500) MeV. For k below k^sB the quarks acquire constituent 
masses Mg ~ 350 MeV due to their Yukawa coupling to the chiral con- 
densate (^. 

• The confinement scale Aqcd — 200 MeV which corresponds to the Lan- 
dau pole in the perturbative evolution of the strong coupling constant a^. 
In our context, this is the scale where possible deviations of the effective 
quark propagator from its classical form and multi-quark interactions 
not included in the meson physics may become very important. 

For scales k in the range k^sB ^ k ^ k^ the most relevant degrees of freedom 
are mesons and quarks. Typically, the dynamics in this range is dominated 



"^One may assume that all other bound states are integrated out. We will comment on this 
issue below. 
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by the strong Yukawa coupling h between quarks and mesons: /i^/(47r) ^ ag. 
One may therefore assume that the dominant QCD effects are included in the 
meson physics and consider a simple model of quarks and mesons only. As 
one evolves to scales below k^sB the Yukawa coupling decreases whereas as 
increases. Of course, getting closer to Aqcd it is no longer justified to neglect 
the QCD effects which go beyond the dynamics of effective meson degrees 
of freedom. On the other hand, the final IR value of the Yukawa coupling 
h is fixed by the typical values of constituent quark masses Mg ~ 350 MeV 
to be /i^/(47r) ~ 4.5. One may therefore speculate that the domination of 
the Yukawa interaction persists down to scales k ~ Mg at which the quarks 
decouple from the evolution of the mesonic degrees of freedom altogether due 
to their mass. Of course, details of the gluonic interactions are expected 
to be crucial for an understanding of quark and gluon confinement. Strong 
interaction effects may dramatically change the momentum dependence of 
the quark n-point functions for k around Aqcd- Yet, as long as one is only 
interested in the dynamics of the mesons one is led to expect that these effects 
are quantitatively no too important. Because of the effective decoupling of the 
quarks and therefore the whole colored sector the details of confinement have 
only little infiuence on the mesonic fiow equations for k <, Aqcd • We conclude 
that there are good prospects that the meson physics can be described by an 
effective action for mesons and quarks for k < k^. The main part of the work 
presented here is concerned with this effective quark meson model. 

In order to obtain this effective action at the compositeness scale fc$ 
from short distance QCD two steps have to be carried out. In a first step 
one computes at the scale kp ~ 1.5 GeV an effective action involving only 
quarks. This step integrates out the gluon degrees of freedom in a "quenched 
approximation" . More precisely, one solves a truncated flow equation for QCD 
with quark and gluon degrees of freedom in presence of an effective infrared 
cutoff kp ~ 1.5 GeV in the quark propagators. The exact flow equation to 
be used for this purpose is obtained by lowering the infrared cutoff Rk for 
the gluons to zero while keeping the one for the quarks flxed. Subsequently, 
the gluons are eliminated by solving the fleld equations for the gluon flelds 
as functionals of the quarks. This will result in a non-trivial momentum 
dependence of the quark propagator and effective non-local four and higher 
quark interactions. Because of the infrared cutoff kp the resulting effective 
action for the quarks resembles closely the one for heavy quarks (at least for 
Euclidean momenta). The dominant effect is the appearance of an effective 
quark potential (similar to the one for the charm quark) which describes the 
effective four-quark interactions. For the effective quark action at kp we only 
retain this four-quark interaction in addition to the two-point function, while 
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neglecting n-point functions involving six and more quarks. 

For typical momenta larger than Aqcd a reliable computation of the ef- 
fective quark action should be possible by usinaJp the quark-giuon flow equa- 
tion a relatively simple truncation. The resultcil far the Fourier transform of 
the potential V{q'^) is shown in fig. H. (See Ref.cil for the relation between 
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Figure 1. Heavy quark potential for different values of a*. The solid line corresponds to 
a* = 1.5 whereas the dashed lines are for «+ = 1 (lower dashed line) and a* = 2 (upper 
dashed line). Also shown are the two loop potential and the Richardson fit as dotted lines. 



the four-quark interaction and the heavy quark potential which involves a 
rescaling in dependence on fcp.) A measure for the truncation uncertainties 
is the parameter a* which corresponds to the value to which an appropri- 
ately defined running strong coupling Us evolves for fc ^ 0. We see that 
these uncertainties affect principally the \ow-q^ region. For high values of 
q^ the potential is very close to the perturbative two-loop potential whereas 
for intermediate q^ relevant for quarkonium spectra it is quite close JiO phe- 
nomenologically acceptable potentials (e.g. the Richardson potential.Ell) The 
inverse quark propagator is found in this computation to remain very well 
approximated by the simple classical momentum dependence ^ . 
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In the second step one has to lower the infrared cutoff in the effective 
non-local quark model in order to extrapolate from kp to fc$. This task can 
be carried out by means of the exact flow equation for quarks only, starting 
at kp with an initial value Tk [V'] as obtained after integrating out the gluons. 
For fermions the trace in (O) has to be replaced by a supeptrace in order to 
account for the minus sign related to Grassmann variables.t2l A first investi- 
gation in this directionEJ has used a truncation with a chirally invariant four 
quark interaction whose most general momentum dependence was retained 



Tk = 










m-''ip)^,+im'^\p)]i,,hip) 



(27r) S{pi +P2-P3- Pi) 



(Pl,P2,P3,P4)| 



i^li-pihsi^iip 



V'fc(P4)V'^(-P3) 



'^fc(P4)75^°(-P3)]}.(44) 



Here i,j run from one to Nc which is the number of quark colors. The in- 
dices a, b label the different light quark flavors and run from 1 to N. The 
matrices m and rh are hermitian and m + ifh'y^ fopais therefore the most 
general quark mass matrix. (Our chiral conventionstj where the hermitian 
part of the mass matrix is multiplied by 75 may be somewhat unusual but 
they are quite convenient for Euclidean calculations.) The ansatz (B) does 
not correspond to the most general chirally invariant four-quark interaction. 
It neglects similar interactions in the p-meson and pomeron channels which 
are also obtained from a Fierz transformation of the heavy quark potential£3 
With V{q^) the heavy quark potential in a Fourier representation, the initial 



value at kr, 



iW 



^r (Pl,P2,P3,P4)^. 



1.5 GeV was taken as [Z^^k 
1, 



2 

■0,fcp 
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-K)) 
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(45) 



This corresponds to an approximation by a one gluon exchange term ~ Ois{kp) 
and a string tension A ~ O.lSGeV^ and is in reasonable agreement with the 
form computed recentlynJ from the solution of flow equations (see figll) . In the 
simplified ansatz (HS) the string tension introduces a second scale in addition 
to kp and it becomes clear that the incorporation of gluon fiuctuations is a 
crucial ingredient_£or the emergence of mesonic bound states. For a more 
precise treatments of the four-quark interaction at the scale fc$ this second 
scale is set by the running of Ug or Aqcd- 
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The evolution equation for the function A^ for k < kp can be derived 



from the fermionic version of (11) and the truncation (|44). Since Aj, depends 



on six independent momentum invariants it is a partial differential equatiaa 
for a function depending on seven variables and has to be solved numericallyJl3 
The "initial value" (EBh corresponds to the i-channel exchange of a "dressed" 
colored gluonic state and it is by far not clear that the evolution of A^ will 
lead at lower scales to a momentum dependence representing the exchange of 
colorless mesonic bound states. Yet, at the compositeness scale 

fc$~630MeV (46) 

one finds an approximate factorization 

Ai.t^bi,P2,P3,P4) = giPi,P2)Gis)g{p3,Pi) + ■■■ (47) 

which indicates the formation of mesonic bound states. Here g{pi,p2) de- 
notes the amputated Bethe-Salpeter wave function and G{s) is the mesonic 
bound state propagator displaying a pole-like structure in the s-channel if it 
is continued to negative s = (pi +^2)^- The dots indicate the part of \^ 
which does not factorize and which will be neglected in the following. In 
the limit where the momentum dependence of g and G is neglected jsje re- 
cover the four-quark interaction of the Nambu-Jona-Lasinio model£3'Ea It is 
therefore not surprising that our description of the dynamics for fc < fc$ will 
parallel certain aspects of other investigations of this model, even though we 
are not bound to the approximations used typically in such studies (large-A^c 
expansion, perturbative renormalization group, etc.). 

It is clear that for scales fc ^ fc$ a description of strong interaction physics 
in terms of quark fields alone would be rather inefficient. Finding physically 
reasonable truncations of the effective average action should be much easier 
once composite fields for the mesons are introduced. The exact renormaliza- 
tion group equation can indeed be supplemented by an exact formalism for 
the introduction of composite field variables or, more generally, a change of 
variables.EJ For our purpose, this amounts in practice to inserting at the scale 
fc$ the identities 

1 ~ I VaA exp \ - tr {(j\ - k\G - m\~ ^"^GX ^ 
X {uA - GKa -mA- go) \ 
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1 ~ / Van exp \ - tr U^ - KJ^G - m\ - O'^'^^'^G) -^ 

X (au - GKh - mn - GO'^^) \ 

into the functional integral which formally defines the quark effective av- 
erage action. Here we have used the shorthand notation A^GB = 
I 7§^^ail)C^'^''i'l)Bbi'l), and Ka,h are sources for the collective fields aA,H 
which correspond in turn to the anti-hcrmitian and hermitian parts of the 
meson field $. They are associated to the fermion bilinear operators 0[ip], 
O^^^lijj] whose Fourier components read 

/74 
J^9{-P,P + q)T{p)MP + q) 

o^'\{q) = -J^g{-P,P + q)T{pHMp + q) ■ (48) 

The choice of g{—p,p + q) as the bound state wave function renormaliza- 
tion and of G(q) as its propagator guarantees that the four-quark interaction 
contained in ( [48[ ) cancels the dominant factorizing part of the QCD-induced 
non-local four-quark interaction Eqs.([44|), (|4^). In addition, one may choose 

m^ab = rhab{0)g-HO,0)Z^,,kAO) (49) 

such that the explicit quark mass term cancels out for q = 0. The remaining 
quark bilinear is ~ m{q)~m{0)Z^^kis,{0)gi—q,q)/[Z^^k,s,iq)g{0,0)]- It vanishes 
for zero momentum and will be neglected in the following. Without loss of 
generality we can take m real and diagonal and m — 0. 

In consequence, we have replaced at the scale fc$ the effective quark ac- 
tion m4) with p7| ) by an effective quark meson action given by 

ffc^Ffc-i /"d4:^tr($tj + /$) 

Tfc^ /"d*a:C/fc($,$t) (50) 

'^''^ > z^M.q)q' tr ['^HqMq)] + z^,k{q)i^aiqh''q,r{q) 



{2^^ 



f d^P T ^ 
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X i^\q) ( ^Ar^'^a.ip) - ^<i>L(-p) ) ^\^-p) 



At the scale k^ the inverse scalar propagator is related to G{q) in ( p7| ) by 

G-\q^)^2ml^+2Z^^kMq^ . (51) 

This fixes the term in C/^^ which is quadratic in $ to be positive, Uk^ = 
rn\ tr$T$ + . . .. The higher order terms in C/^^ cannot be determined in 
the approximation (^4|) since they correspond to terms involving six or more 
quark fields. The initial value of the Yukawa coupling corresponds to the 
"quark wave function in the meson" in (|47|), i.e. 

hk^{-q,q-p) = 9{-q,q-p) (52) 

which can be normalized with /ife^(0,0) = 5(0,0) = 1. We observe that the 
explicit chiral symmetry breaking from non-vanishing current quark masses 
appears now in the form of a meson source term with 

J = 2mfc^Z^.fc^(0)g"^(0,0) {mab + irhab) = 2Z^.fc^m^_j,diag(m„, Wd, m^) . 

(53) 
This induces a non-vanishing ($) and an effective quark mass Mq through 
the Yukawa coupling. We note that the current quark mass mq and the con- 
stituent quark mass Mq ~ hk ($) are identical at the scale fc$. (By solving 
the field equation for $ as a functional of ip, tp (with Uk — m|tr$t(j)) one 
recovers from ( |50|) the effective quark action (^^. For a generalization be- 
yond thjC-japproximation of a four-quark interaction or a quadratic potential 
see Ref.Ea.) Spontaneous chiral symmetry breaking can be described in this 
language by a non-vanishing ($) in the limit j — > 0. Because of spontaneous 
chiral symmetry breaking the constituent quark mass Mq can differ from zero 
even for mq = 0. It is a nice feature of our formalism that it provides for a 
unified description of the concepts of the current and the constituent quark 
masses. As long as the effective average potential Uk has a unique minimum 
at $ = there is simply no difference between the two. The running of the 
current quark mass in the pure quark model should be equivalent in the quark 
meson language to the running of hk (<&) (k). (A verification of this property 
would actually provide a good check for the truncation errors.) Nevertheless, 
the formalism is now adapted to account for the quark mass contribution from 
chiral symmetry breaking since the absolute minimum of Uk niay be far from 
the perturbative one for small k. We also note that Tk in (pO) is chirally 
invariant. The explicit chiral symmetry breaking in Tk appears only in the 
form of a linear source term which is independent of k and does not affect the 
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flow of Ffe. The flow equation for Tk therefore respects chiral symmetry even 
in the presence of quark masses. This leads to a considerable simplification. 
At the scale fc$ the propagator G and the wave function g{—q, q — p) 
should be optimized for a most complete elimination of terms quartic in the 
quark fields. In the present context we will, however, neglect-±hc momentum 
dependence of Z^^k, Z^^k and h^. The mass toj;^ was found inCJ for the simple 
truncation (H) with Z-^^^ = 1, m = to = to be m^^ ~ 120 MeV. In view of 
the possible large truncation errors we will take this only as an order of mag- 
nitude estimate. Below we will consider the range mt^ — (45 — 120) MeV for 
which chiral symmetry breaking can be obtained in a two fiavor model. Fur- 
thermore, we will assume, as usually done in large-A'c computations within 
the NJL-model, that Z^^k,^ = Z^^k^(q = 0) <C 1. The quark wave function 
renormalization Z^^k = Z^^k{q = 0) is set to one at the scale k^ for conve- 
nience. For k < k^ we will therefore study an effective action for quarks and 
mesons in the truncation 

Ffc = Jd^xl Z^,fe^„z^ V"^ + ^$,fctr [d,,<i>^d^^] + Uk{^, $+) 

+ h.r (^*- - ^^(<i>^)a.) A (54) 

with compositeness conditions 
Uk^ ($, $t) = ml^ tr $^$ - -Vk^ (det $ -I- det $t) 

+ ^\m (tr$^$)%^— A2,fc^trf$^$-— tr*^* 
^2 - ^ -- (45MeV)2- (120MeV)2 



20(0) 

hk^ — Z^^kg, — 1 

^*,fc. « 1 . (55) 

As a consequence, the initial value of the renormalized Yukawa coupling which 
is given by hk^ = hkis,ZZ,\ Z^ ^ is large! Note that we have included in the 
potential an explicit /7a (1) breaking term H ~ F^ which mimics the effect of 
the chiral anomaly of QCD to leading order in an expansion of the effective 
potential in powers of <&. Because of the infrared stability of the evolution of F^ 

■'^The anomaly term in the fermionic effective average action has been computed in.P^ 
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which will be discussed below the precise form of the potential, i.e. the values 
of the quartic couplings A^./t^ and so on, will turn out to be unimportant. 

We have refrained here for simplicity from considering four quark opera- 
tors with vector and pseudo-vector spin structure. Their inclusion is straight- 
forward and would lead to vector and pseudo-vector mesons in the effective 
action (|3). We will concentrate first on two flavors and consider only the 
two limiting cases F^ = and F^ — + cx3. We also omit first the explicit quark 
masses and study the chiral limit j = 0. Because of the positive mass term 
to| one has at the scale fc$ a vanishing expectation value ($) = (for j = 0). 
There is no spontaneous chiral symmetry breaking at the compositeness scale. 
This means that the mesonic bound states at fc$ and somewhat below are not 
directly connected to chiral symmetry breaking. 

The question remains how chiral symmetry is broken. We will try to 
answer it by following the evolution of the effective potential Uk from fc$ to 
lower scales using the exact renormalization group method outlined above with 
the compositeness conditions ( |55| ) defining the initiaLvalues. In this context it 
is important that the formalism for composite field^ also induces an infrared 
cutoff in the meson propagator. The flow equations are therefore exactly of the 
form (|l|) (except for the supertrace), with quarks and mesons treated on an 
equal footing. In fact, one would expect that the large renormalized Yukawa 
coupling will rapid|ly, drive the scalar mass term to negative values as the IR 
cutoff k is lowered. t2l This will then finally lead to a potential minimum away 
from the origin at some scale fc^sB < k^ such that ($) ^ 0. The ultimate goal 
of such a procedure, besides from establishing the onset of chiral symmetry 
breaking, would be to extract phenomenological quantities, like f-^ or meson 
masses, which can be computed in a straightforward manner from Tk in the 
IR hmit fc ^ 0. 

At first sight, a reliable computation of Ffe^o seems a very difficult task. 
Without a truncation Tk is described by an infinite number of parameters 
(couplings, wave function renormalizations, etc.) as can be seen if Tk is ex- 
panded in powers of fields and derivatives. For instance, the pseudoscalar 
and scalar meson masses are obtained as the poles of the exact propaga- 

(2) 

tor, lim/c^o F^ ((7)|$=($), which receives formally contributions from terms in 
Ffe with arbitrarily high powers of derivatives and the expectation value ctq. 
Realistic nonperturbative truncations of F^ which reduce the problem to a 
manageable size are crucial. We will argue in the following that there may be 
a twofold solution to this problem: 

• Due to an IR fixed point structure of the flow equations in the symmetric 
regime, i.e. for k^sB < k < k^, the values of many parameters of F^ 
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for fc ^ will be approximately independent of their initial values at 
the compositeness scale A:$. For small enough Z$.fc^ only a few relevant 
parameters (m^ , Ffe^) need to be computed accurately from QCD. They 
can alternatively be determined from phenomenology 

• One can show that physical observables like meson masses, decay con- 
stants, etc., can be [Expanded in powers of the quark masses within the 
linear meson model.E3'E3'E3jrhis is similar to the way it is usually done in 
chiral perturbation theoryEa To a given finite order of this expansion only 
a finite number of terms of a simultaneous expansion of Ffc in powers of 
derivatives and $ are required if the expansion point is chosen properly. 

In combination, these two results open the possibility for a perhaps unex- 
pected degree of predictive power within the linear meson model. 

We wish to stress, though, that a perturbative treatment of the model 
at hand, e.g., using perturbative RG techniques, cannot be expected to yield 
reliable results. The renormalized Yukawa coupling is expected to be large at 
the scale kg, . Even the IR value of h is still relatively big 

lini h ^ —A ^ 7.5 (56) 

and h increases with increasing k. The dynamics of the linear meson model 
is therefore clearly nonperturbative for all scales k <k^. 

6 Flow equations for the linear quark meson model 

We will next turn to the ERGE analysis of the linear meson model which 
was introduced in the last section^. In a first approach we will attack the 
problem at hand by truncating F^ in such a way that it contains all perturba- 
tively relevant and marginal operators, i.e. those with canonical dimensions. 
This has the advantage that the flow of a small number of couplings permits 
quantitative insight into the relevant mechanisms, e.g., of chiral symmetry 
breaking. More quantitative precision will be obtained once we generalize our 
truncation (see below) for the 0(4)-model by allowing for the most general 
effective average potential Uk- The effective potential Uk is a function of only 
four SUl{N) X SUr{N) invariants for A^ = 3: 

p = tr$t$ 

''For a study of chiijftl symmetry breaking in QED using related exact renormalization group 
techniques see Ref.Ej. 



Jaro: submitted to World Scientific on December 18, 1998 27 



-. = ^tr(<i>t$_lp' 



1 ^^ 



T3 = tr ( $^$ 

C = det $ + det <^^ . (57) 

The invariant T3 is only independent for A^ > 3. For A^ = 2 it can be elimi- 
nated by a suitable combination of T2 and p. The additional C/a(1) breaking 
invariant w = i(det $ — det $^) is CV violating and may therefore appear only 
quadratically in Uk- It is straightforward to see that w^ is expressible in terms 
of the invariants ( pTj) . We may expand Uk as a function of these invariants 
around its minimum, i.e. p = po = -^^O' C = Co = '^^q E^nd T2 = T3 = where 

ao^^tr($). (58) 

The expansion coefficients are the fc-dependent couplings of the model. In 
our first version we only keep couplings of canonical dimension dc < 4. This 
yields in the chirally symmetric regime, i.e., for fc^sB < /c < fc$ where ao = 

Uk = rnlp + -Xi,kP^ H r — X2,kT2 - -VkS, (59) 

whereas in the SSB regime for k < /c^sb we have 

Uk = \\iM [P - Nalk] ' + ^^ A2,fcr2 + \vk [<,7 V - C] • (60) 

Before continuing to compute the nonperturbative beta functions for these 
couplings it is worthwhile to pause here and emphasize that naively (pertur- 
batively) irrelevant operators can by no means always be neglected. The most 
prominent example for this is QCD itself. It is the very assumption-of our 
treatment of chiral symmetry breaking (substantiated by the results oiE2l) that 
the momentum dependence of the coupling constants of some six-dimensional 
quark operators (i/'V')^ develop poles in the s-channel indicating the formation 
of mesonic bound states. On the other hand, it is quite natural to assume that 
(j)6 or $^ operators are not really necessary to understand the properties of 
the potential in a neighborhood around its minimum. Yet, truncating higher 
dimensional operators does not imply the assumption that the corresponding 
coupling constants are small. In fact, this could only be expected as long 
as the relevant and marginal couplings are small as well. What is required, 
though, is that their influence on the evolution of those couplings kept in the 
truncation, for instance, the set of equations ( |62| ) below, is small. A com- 
parison with the results of the later sections for the full potential Uk (i.e.. 
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including arbitrarily many couplings in the formal expansion) will provide a 
good check for the validity of this assumption. In this context it is perhaps 
also interesting to note that the truncation (M) includes the known one-loop 
beta functions of a small coupling expansion as well as the leading order re- 
sult of the large- A^c expansion of the Ul{N) x Ur{N) modelS This should 
provide at least some minimal control over this truncation, even though we 
believe that our results are significantly more accurate. 

Inserting the truncation Eqs.(|4|), (p9[), ( |60[ ) into ( pi] ) reduces this func- 
tional differential equation for infinitely many variables to a finite set of ordi- 
nary differential equations. This yields, in particular, the beta functions for 
the couplings A]_fe., A2,fe, f^fc and rrTJ. or CTo.fc- Details of the calculation can 
be found in Ref.E3. We will refrain here from presenting the full set of flow 
equations but rather illustrate the main results with a few examples. Defining 
dimensionless renormalized VEV and coupling constants 






Ai — XihZ, 



i,^'^<^.k ' 



i = 1,2 



V = VkZ^ I k 



N-4 



(61) 



one finds, e.g., for the spontaneous symmetry breaking (SSB) regime and 
77,. = n 



Vk = 



dt 



n = —(2 + ??$)« + 



167r2 



N^l\{Q-m) + il\{2\iK-r^^) 



+ {N^ - 1) 



-^ 



/f(A2K;r,<,)-4iV,^/p(l/i2«;;^^) 



^Ai = 27y<E.Ai + :^ \ N^XlliiO; m) + 9A?Z|(2Aik; 77$) 



dt 



167r2 



+ {N' - 1) [Ai + X,f II{\2k; ry$) - i^h¥/^\^h^n; rj^) \ (62) 



1 iV2 

167r^ 4 •^ 



d 
dt' 



9 

4^ 



-N'^XJll^iO, X2H; m) + 3[A2 + 4Ai]A2/ti(2AiK, A2K; m) 
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^/i2 ^ [d - 4 + 2r,^ + 77$] h^ - ^Vdh^ |iv2z(^,^)''(l«;/i^ 0; ij^,m) 

- '1,1 ']v'* '2kAi; 7/^,77$) 

Here 77$ — —4; InZ^^k, iM,~ ~li ^^^^,k are the meson and quark anomalous 
dimensions, respectively^ The symbols if^, Ini^n^ ^^"^ '" denote bosonic apd 
fermionic mass threshold functions, respectively, which are defined in Ref.H. 
They describe the decoupling of massive modes and provide an important 
nonperturbative ingredient. For instance, the bosonic threshold functions 

"^ '^ 4 "^ J {2T:YZ^,k dt [P(92) + p^,]»+i ^ ^ 

involve the inverse average propagator P(g^) = q^ + Z^\Rk{q^) where the 
infrared cutoff is manifest. These functions decrease ~ ui~'-"+^-' for w ^ 1. 
Since typically w = Al'^/k'^ with M a mass of the model, the main effect 
of the threshold functions is to cut off fluctuations of particles with masses 
M^ ^ k"^. Once the scale k is changed below a certain mass threshold, the 
corresponding particle no longer contributes to the evolution and decouples 
smoothly. 

Within our truncation the beta functions ( plj) for the dimensionless cou- 
plings look almost the same as in one-loop perturbation theory. There are, 
however, two major new ingredients which are crucial for our approach: First, 
there is a new equation for the running of the mass term in the symmet- 
ric regime or for the running of the potential minimum in the regime with 
spontaneous symmetry breaking. This equation is related to the quadratic 
divergence of the mass, term in perturbation theory and does not appear in 
the Callan-SymanziliO or Coleman- Weinbergp^ treatment of the renormal- 
ization group. Obviously, these equations are the key for a study of the onset 
of spontaneous chiral symmetry breaking as k is lowered from fc$ to zero. 
Second, the most important nonperturbative ingredient in the flow equations 
for the dimensionless Yukawa and scalar couplings is the appearance of effec- 
tive mass threshold functions like (p3) which account for the decoupling of 
modes with masses larger than k. Their form is different for the symmetric 
regime (massless fermions, massive scalars) or the regime with spontaneous 
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symmetry breaking (massive fermions, massless Goldstone bosons). Without 
the inclusion of the threshold effects the running of the couplings would never 
stop and no sensible limit /c — > could be obtained because of unphysical 
infrared divergences. The threshold functions are not arbitrary but have to 
be computed carefully. The mass terms appearing in these functions involve 
the dimcnsionless couplings. Expanding the threshold functions in powers of 
the mass terms (or the dimensionless couplings) makes their non-perturbative 
content immediately visible. It is these threshold functions which will make it 
possible below to use one-loop type formulae for the necessarily nonperturba- 
tive computation of the critical behavior of the (effectively) three-dimensional 
0(4)-symmetric scalar model. 

7 The chiral anomaly and the 0(4) model 

We have seen how the mass threshold functions in the flow equations describe 
the decoupling of heavy modes from the evolution of Tk as the IR cutoff k is 
lowered. In the chiral limit with two massless quark flavors (N — 2) the pions 
are the massless Goldstone bosons. Below, once we will consider the flow of the 
full effective average potential without resorting to a quartic truncation, we 
will also consider the more general case of non- vanishing current quark masses. 
For the time being, the effect of the physical pion mass of ni.^ — 140 MeV, 
or equivalently of the two small but non-vanishing current quark masses, can 
easily be mimicked by stopping the flow of T^ at k = tti^ by hand. This 
situation changes significantly once the strange quark is included. Now the 
ry and the four K mesons appear as additional massless Goldstone modes 
in the spectrum. They would artificially drive the running of Tk at scales 
rriTr ^ k ^ 500 MeV where they should already be decoupled because of their 
physical masses. It is therefore advisable to focus on the two flavor case N = 2 
as long as the chiral limit of vanishing current quark masses is considered. 

It is straightforward to obtain an estimate of the (renormalized) coupling 
ly which parameterizes the explicit C/yi(l) breaking due to the chiral anomaly. 
From ( |60| ) we find 

ml^^iya^-^^lGeY (64) 

which translates for A^ = 2 for fc ^ into 

i^{k ^ 0) ~ 1 GeV . (65) 

This suggests that j/ — > cx) can be considered as a realistic limit. An important 
simplification occurs for iV = 2 and z^ — > oo, related to the fact that for A^ = 2 
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the chiral group SUl{'2) x SUniT) is (locally) isomorphic to 0(4). Thus, the 
complex (2, 2) representation $ of SUl{'2) x SUii{2) may be decomposed into 
two vector representations, (cr, tt'^) and {r]',a'^) of 0(4): 

* = ^ (^ - w/) + ^ (a' + *^') ^k ■ (66) 

For 1/ ^ oo the masses of the rf and the a*^ diverge and these particles 
decouple. We are thenJsft with the original 0(4) symmetric linear cr-model 
of Gell-Mann and LevyE3 coimled to quarks. The flow equations of this model 
have been derived previouslyHQ for the truncation of the effective action used 
here. For mere comparison we also consider the opposite limit i^ — > 0. Here 
the rj' meson becomes an additional Goldstone boson in the chiral limit which 
suffers from the same problem as the K and the ry in the case TV = 3. Hence, 
we may compare the results for two different approximate limits of the effects 
of the chiral anomaly: 

• the 0(4) model corresponding to A^ = 2 and i^ ^ oo 

• the C/l(2) X Ufi{2) model corresponding to A^ = 2 and v = Q. 

For the reasons given above we expect the first situation to be closer to reality. 
In this case we may imagine that the fluctuations of the kaons, ry, rj' and 
the scalar mesons (as well as vector and pseudovector mesons) have been 
integrated out in order to obtain the initial values of F^.^ — in close analogy 
to the integration of the gluons for the effective quark action F^; [ip] discussed 
above. We will keep the initial values of the couplings rn\ , Ai^^^ and Z^^kis, as 
free parameters. Our results should be quantitatively accurate to the extent 
to which the local polynomial truncation is a good approximation. 

8 Infrared stability 

Dz and the corresponding set of flow equations for the symmetric regime con- 
stitute a coupled system of ordinary differential equations which can be inte- 
grated numerically. Similar equations can be computed for the 0(4) model 
where N — 2 and the coupling A2 is absent. We have neglected in a first step 
the dependence of all threshold functions appearing in the flow equations on 
the anomalous dimensions. This dependence will be taken into account below 
once we abandon the quartic potential approximation. The most important 
result is that chiral symmetry breaking indeed occurs for a wide range of ini- 
tial values of the parameters including the presumably realistic case of large 
renormalized Yukawa coupling and a bare mass rrlfe^ of order 100 MeV. A 
typical evolution of the renormalized meson mass m with k is plotted in fig. y. 
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Figure 2. Evolution of the rcnormalizcd mass m in the symmetric regime (dashed line) and 
the vacuum expectation value ctq = Z',o-qi^ of the scalar field in the SSB regime (solid 
line) as functions of k for the Ul{2) x Uji{2) model. Initial values at A; = fc$ are A; / = 0, 
h] = 300 and nik^ = 63 MeV. 



Driven by the strong Yukawa coupling, m decreases rapidly and goes through 
zero at a scale fc^sB not far below fc$. Here the system enters the SSB regime 
and a non-vanishing (renormalized) VEV (Tq for the meson field $ develops. 
The evolution of ctq with k turns out to be reasonably stable already before 
scales k ~ TOtt where the evolution is stopped. We take this result as an 
indication that our truncation of the effective action Fj. leads at least qualita- 
tively to a satisfactory description of chiral symmetry breaking. The reason 
for the relative stability of the IR behavior of the VEV (and all other cou- 
plings) is that the quarks acquire a constituent mass Mq — ha^ ~ 350 MeV 
in the SSB regime. As a consequence they decouple once k becomes smaller 
than Mq and the evolution is then dominantly driven by the massless Gold- 
stone bosons. This is also important in view of potential confinement effects 
expected to become important for the quark n-point functions for k around 
Aqcd — 200 MeV. Since confinement is not explicitly included in our trun- 
cation of Ffe, one might be worried that such effects could spoil our results 
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completely. Yet, as discussed in some more detail above, only the colored 
quarks should feel confinement and they are no longer important for the evo- 
lution of the meson couplings for k around 200 MeV. One might therefore 
hope that a precise treatment of confinement is not crucial for this approach 
to chiral symmetry breaking. 

Most importantly, one finds that the system of flow equations exhibits a 
partial IR fixed point in the symmetric phase. As already pointed out one 
expects Zi^^kis, to be rather small. In turn, one may assume that, at least for the 
initial range of running in the symmetric regime the mass parameter m^ ~ 
Z^ 1^ is large. This means, in particular, that all threshold functions with 



arguments 



may be neglected in this regime. As a consequence, the fiow 



equations simplify considerably. We find, for instance, for the C/i(2) x Ur{2) 
model 



(67) 





d Xi ^ 
Jtl?~ 




f^^- 
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dt 


d X2 ^ 
Jtl^~ 
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'^ 7,2 ^ 
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m 



87r2 
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This system possesses an attractive IR fixed point for the quartic scalar self 
interactions 



A -^A - 2 
Ai. - -\2. - - 



(68) 



Furthermore it is exactly soluble.!^ Because of the strong Yukawa coupling 
the quartic couplings Ai and A2 generally approach their fixed point values 
rapidly, long before the systems enters the broken phase (to — > 0) and the 
approximation of large to breaks down. In addition, for large initial values /ij 
the Yukawa coupling at the scale k^sB where to vanishes (or becomes small) 
only depends on the initial value mfe^. Hence, the system is approximately 
independent in the IR upon the initial values of Ai, A2 and /i^, the only 
"relevant" parameter being rfik^ . (Once quark masses and a proper treatment 
of the chiral anomaly are included for iV = 3 one expects that ruq and u are 
additional relevant parameters. Their values may be fixed by using the masses 
of TT, K and 77' as phenomenological input.) In other words, the effective action 
looses almost all its "memory" in the far IR of where in the UV it came from. 
This feature of the flow equations leads to a perhaps surprising degree of 
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Figure 3. The pion decay constant /,r as a function of m| /fc|, for fccj = 630 MeV and 
initial values A;,/ = and /i| = 300 (solid line) as well as /i| = lO** (dashed line). 



predictive power. In addition, also the dependence of /jr — 2ao on 771^^ is not 
very strong for a large range in rfik^ , as shown in fig. || 

The relevant parameter Wife^ can be fixed by using the constituent quark 
mass Mq = hao ~ 350 MeV as a phenomenological input. One obtains for the 
0{4)-model 



ttt2 



J^2 



0.02. 



(69) 



The resulting value for the decay constant is 



/^ = (91 - 100) MeV (70) 

for hj — 10* — 300. It is striking that this comes close to the real value 
/tt — 92.4 MeV but we expect that the uncertainty in the determination of 
the compositeness scale fc$ and the truncation errors exceed the influence of 
the variation of h^. We have furthermore used this result for an estimate of 
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the chiral condensate: 

(V^V) = -^L /-%l=o^ - -(195MeV)3 (71) 

where the factor A ~ 1.7 accounts for the change of the normahzation scale of 
{^V) from fc$ to the commonly used value 1 GeV. Our value is in reasonable 
agreement with results from sum rules.E3 This is non-trivial since not only 
mfc^ and /^ enter but also the IR value Z,^k=o- Integrating ( pTj ) for 77$ one 
finds 

^*,fc = ^*,fc„, + 7—7 In -r ■ (72) 

Stt^ k 

Thus Z^js will indeed be practically independent of its initial value Z^^k^ 
already after some running as long as 2'$^^^ is small compared to 0.01. 

The alert reader may have noticed that the beta functions (|6^) corre- 
spond exactly to those obtained in the one-quark-loop approximation or, in 
other words, to the leading order in the large-iVc expansion for the Nambu- 
Jona-Lasinio model.EHl The fixed point ( pq ) is then nothing but the large-Ai'c 
boundary condition on the evolution of Ai and A2 in this model. Yet, we wish 
to stress that nowhere we have made the assumption that Nc is a large number. 
On the contrary, the physical value Nc = 3 suggests that the large-TV^ expan- 
sion should a priori only be trusted on a quantitatively rather crude level. The 
reason why we expect ( pT\ ) to nevertheless give rather reliable results is based 
on the fact that for small Z^^k,;, all (renormalized) meson masses are much 
larger than the scale k for the initial part of the running. This implies that the 
mesons are effectively decoupled and their contribution to the beta functions 
is negligible leading to the one-quark-loop approximation. Yet, already after 
some relatively short period of running the renormalized meson masses ap- 
proach zero and our approximation of neglecting mesonic threshold functions 
breaks down. Hence, the one-quark-loop approximation is reasonable only 
for scales close to fc* but is bound to become inaccurate around fc^sB and in 
the SSB regime. What is important in our context is not the numerical value 
of the partial fixed points (^8|) but rather their mere existence and the pres- 
ence of a large coupling h"^ driving the A^ fast towards them. This is enough 
for the IR values of all couplings to become almost independent of the initial 
values Xij. Similar features of IR stability are expected if the truncation is 
enlarged, for instance, to a more general form of the effective potential C/fc<^ 
as will be discussed in the next section. 
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9 Flow equation for the scalar potential 

In this and the remaining sections we consider the 0(4)-syninietric quark me- 
son rondel without truncating its effective average potential to a polynomial 
form.Ea A comparison with the results of the last sections will give us a feel- 
ing for the size of the errors induced by the quartic truncation used so far. 
Furthermore, such an approach is well suited for a study of the chiral phase 
transition close to which the form of the potential deviates substantially from 
a polynomial. 

It is convenient to work with dimensionless and renormalized variables 
therefore eliminating all explicit fc-depcndence. With {t — \n{k/k<j,)) 

u{t, p) = k-^'Ukip) , P = Z^,kk^'^P (73) 

and using (p4) as a first truncation of the effective average action Tk one 
obtains the flow equation in arbitrary dimensions d 

d 

—-U ~ —du + {d~ 2 + 77$) pu' 
at 

+ 2vd Ult,{u; m) + li{u + 2fm"- m) - 2i+^N/p\]^ph'-T^^) 

^(74) 

Here v^^ = 2''+^7r''/^r((i/2) and primes denote derivatives with respect to 
p. We will always use in the following for the number of quark colors N^ = 
3. dTJ) is a partial differential equation for the effective potential u{t, p) 
which has to be supplemented by the flow equation for the Yukawa coupling 
and expressions for the anomalous dimensions 77*, 77^ . The definition of the 

threshold functions l^^ In can be found in Ref.S 

The dimensionless renormalized expectation value k = 2fc^^'^Z$^fcCTQ ;., 
with CTo.fc the fc-dependent VEV of $, may be computed for each k directly 
from the condition 

u'{t, k) = -J=k-"-^ Z^]i^ = eg (75) 



/2k 

whereE3 

J = 2m2^,fem (76) 

and m = {mu+md)/2 denotes the average light current quark mass normalized 
at fc$. Note that k = in the symmetric regime for vanishing source term. ( |75|) 
allows us to follow the flow of k according to 
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with A = u"{t—^). We define the Yukawa coiiphng for p ~ k and its flow 
equation readsL^ 

— h^ = (d - 4 + 2r/^ + ?7$) K^ - 2vdh'^l 3/^ ^ ii^h'^K, eg-, ri^,m) 

Similarly, the the scalar and quark anomalous dimensions are inferred from 
m = -— lnZ$,fc = 4— •! 4KA^m2_2(es,eg + 2Ak;?7$) 



(79) 



Vip = --irln^v.fc = 2— ft < 3m^2 (^" n,<^g;'n^,'n'S>J 

{FB)dA ,2 , r,\ ^ I 

+ "'-b (2"' '«, £3 + 2Ak; 77^,77$) ^ , 

which is a linear set of equations for the anomalous dimensions. The defini- 
tions of the threshold functions lni,n2 , f^n^ n^i ^4 ^nd mni.ni Siie again 
specified in Ref E^ . 

The flow equations (^J), ( [77| ) — (|79|), constitute a coupled system 
of ordinaryj— fLnd partial differential equations which can be integrated 
numericallyO' l3 Here we take the effective current quark mass dependence of 
h, Z^^k and Z^,k into account by stopping the evolution according to Eqs.dTq), 
([79|), evaluated for the chiral limit, below the pion mass TOtt- 

Similarly to the case of the quartic truncation of the effective average 
potential described in the preceding section one finds for d = 4 that chiral 
symmetry breaking indeed occurs for a wide range of initial values of the 
parameters. These include the presumably realistic case of large renormalizcd 
Yukawa coupling and a bare mass mfc^ of order 100 MeV. 

Most importantly, the approximate partial IR fixed point behavior en- 
countered for the quartic potential approximation before, carries over to. the 
truncation of Tk which maintains the full effective average potentials To 
see this explicitly we study the flow equations (fz^), ( ^T\ ) — ( [79| ) subject to 
the condition Z^^k^ <C 1. For the relevant range of p both u'(t,p) and 
u'{t,p) -\- 2pu"{t,j5) are then much larger than ph?(t) and we may there- 
fore neglect in the flow equations all scalar contributions with threshold 
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functions involving these large masses. This yields the simplified equations 
(d = 4, v^^ = 327r2) 

Again, this approximation is only valid for the initial range of running below 
fc$ before the (dimensionless) renormalized scalar mass squared u'{t,p = 0) 
approaches zero near the chiral symmetry breaking scale. The system ( pO|) is 
exactly soluble. We find 



1 - i^hjt ' 
u{t,p) = e-'^ujie^'p^) - ^ iyre-'rip\lh\t)pe'n ■ 



(81) 



(The integration over r on the right hand side of the solution for u can be 
carried out by first exchanging it with the one over momentum implicit in 
the definition of the threshold function Zq .) Here u/(/5) = u{0,p) denotes 
the effective average potential at the compositeness scale and hj is the initial 
value of ft.^ at fc^, i.e. for t = 0. For simplicity we will use an expansion of 
the initial value effective potential ui (p) in powers of p around /5 = 



u 



(n 



r(o),„ 



M~P)-T.-^P- (82) 

n=0 

even though this is not essential for the forthcoming reasoning. Expanding 

,(F)4 



also /q in ( pl[ ) in powers of its argument one finds for n > 2 



h^^it) /if '^ ^2 2-+2(n-2) " ^ ' 



1-e 



2(n-2)t 



(83) 

For decreasing i ^ — oo the initial values u^ become rapidly unimportant 
and m(")//i^" approaches a fixed point. For n — 2, i.e., for the quartic couphng, 
one finds 



u(')(t,0) 



1 - . nA.. (84) 
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leading to the fixed point value (u^^V^^)* = 1 already encountered in (p8|). As 
a consequence of this fixed point behavior the system looses all its "memory" 
on the initial values Uj- at the compositeness scale kii>\ This typically 
happens before the approximation u'(t, p), u'(t, p)+2pu"{t, p) ^ ph'^{t) breaks 
down and the solution ( pl| ) becomes invalid. Furthermore, the attraction to 
partial infrared fixed points continues also for the range of k where the scalar 
fluctuations cannot be neglected anymore. As for the quartic truncation the 
initial value for the bare dimensionless mass parameter 



hj k^ 

is never negligible. (In fact, using the values for to|^ and fc$ computed 
previousl}^ for a pure quark effective action as described above one obtains 
rn\^/k% ~ 0.036.) For large hi (and dropping the constant piece ui{Q)) the 
solution (^l]) therefore behaves with growing \t\ as 



^*W--^i, ^^(t) = i 



87r2 
h'it) ^ -— , (86) 



Nrt 



In other words, for /i/ — > cx) the IR behavior of the linear quark meson model 
will depend (in addition to the value of the compositeness scale fc$ and the 
quark mass m) only on one parameter, m| . We have numerically verified 

(2'] 

this feature by starting with different values for Uj (0). Indeed, the differ- 
ences in the physical observables were found to be small. This IR stability of 
the flow equations leads to a perhaps surprising degree of predictive power: 
Not only the scalar wave function rcnormalization but even the full effective 
potential U{p) — limfc_^o Uk{p) is (approximately) fixed for Z^^/t^ <C 1 once 
frikiS! is known! For definiteness we will perform our numerical analysis of the 
full system of flow equations (jT^), (77) — ( |79| ) with the idealized initial value 



ui{p) = u'j{0)p in the limit hj —i- oo. It should be stressed, though, that 
deviations from this idealization will lead only to small numerical deviations 



. inear quark meson model as long as the condition 
hi ^ 15. 

With this knowledge at hand we may now flx the remaining three param- 
eters of our model, kg,, rnl. and rh by using f.^ — 92.4 MeV, the pion mass 



in the IR behavior of th 
■^$,fe<t ^ 1 holds, say forE^ 
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Table 1. The table shows the dependence on the constituent quark mass Mq of the input 



parameters fc$, mf, /k^ 

input used here besides Mq is /,r 



92.4 MeV, m^ 



135 MeV. 



The phenomenological 
The first line corresponds 



demonstrate the insensitivity of our results with respect to the precise values of these 
parameters. 



McV 




fc.,. 

McV 




McV 


m(lGcV) 
McV 


(•0V)(lGeV) 
McV^ 


flO) 

McV 


303 


1 


618 


0.0265 


14.7 


11.4 


-(186)^ 


80.8 


300 





602 


0.026 


15.8 


12.0 


-(183)^ 


80.2 


310 





585 


0.025 


16.9 


12.5 


-(180)^ 


80.5 


339 





552 


0.0225 


19.5 


13.7 


-(174)^ 


81.4 



M^ — 135 MeV and the constituent quark mass Mq as phenomenological in- 
put. This approach differs from that of the preceding sections where we took 
an earlier determination of /c$ as input for the computation of /tt- It will be 
better suited for precision estimates of the high temperature behavior in the 
following sections. Because of the uncertainty regarding the precise value of 
Mq we give in tab. |l] the results for several values of Mq. The first line of 
tab. |l] corresponds to the choice of Mq and Xj = u'j{k) which we will use 
for the forthcoming analysis of the model at finite temperature. As argued 
analytically above the dependence on the value of A/ is weak for large enough 
hi as demonstrated numerically by the second line. Moreover, we notice that 
our results, and in particular the value of j, are rather insensitive with respect 
to the precise value of Mq. It is remarkable,that the values for fc$ and m^^ 
are not very different from those computedc2l from four -quark interactions as 
described above. As compared to the analysis of the preceding sections the 
present truncation of T^ is of a higher level of accuracy: We now consider 
an arbitrary form of the effective average potential instead of a polynomial 
approximation and we have included the pieces in the threshold functions 
which are proportional to the anomalous dimensions. It is encouraging that 
the results are rather robust with respect to these improvements. 



Once the parameters fc$, 



and 771 are fixed there are a number of 



"predictions" of the linear meson model which can be compared with the 
results obtained by other methods or direct experimental observation. First 
of all one may compute the value of jfi at a scale of 1 GeV which is suitable for 
compiaxison with results obtained from chiral perturbation theoryO and sum 
rules O For this purpose one has to account for the running of this quantity 
with the normalization scale from fc$ as given in tab. ^ to the commonly 
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used value of 1 GeV: m(l GeV) = A~^Th{k<si). A reasonable estimate of t 
factor A is obtained from the three loop running of rh in the MS scheme 
For Mq ~ 300 MeV corresponding to the first two lines in tab. its value is 
A ~ 1.3. The results for miLGeV) are in acceptable agreement with recent 
results from other methodsEirEJ even though they tend to be somewhat larger. 
Closely related to this is the value of the chiral condensate 

(V^V) (1 GeV) ^ -Ami [/-^*i=o - 2™] • (87) 

These results are quite non-trivial since not only f-^ and rn^ enter but also 
the computed IR value Z,^^k=o- We emphasize in this context that there 
may be substantial corrections both in the extrapolation from fc$ to 1 GeV 
and because of the neglected influence of the strange quark which may be 
important near A:$. These uncertainties have only little effect on the physics 
at lower scales as relevant for our analysis of the temperature effects. Only 
the value of j which is fixed by ?n,r enters here. 

A further more qualitative test concerns the mass of the sigma resonance 
or radial mode in the limit k —f whose renormalized mass squared is given 

by 



From our numerical analysis we obtain A^^o — 20 which translates into m^ ~ 
430 MeV. One should note, though, that this result is presumably not very 
accurate as we have employed in this work the approximation of using the 
Goldstone boson wave function renormalization constant also for the radial 
mode. Furthermore, the explicit chiral symmetry breaking contribution to 
rn^ is certainly underestimated as long as the strange quark is neglected. In 
any case, we observe that the sigma meson is significantly heavier than the 
pions. This is a crucial consistency check for the linear quark meson model. 
A low sigma mass wmild be in conflict with the numerous successes of chiral 
perturbation theoryo which requires the decoupling of all modes other than 
the Goldstone bosons in the IR-limit of QCD. The decoupling of the sigma 
meson is, of course, equivalent to the limit A — + cxd which formally describes the 
transition from the linear to the non-linear sigma model and which appears 
to be reasonably well realized by the large IR~values of A obtained in our 
analysis. We also note that the issue of the sigma mass is closely connected 
to the value of /i , the value of f-^ in the chiral limit rh — Q also given in 
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Figure 4. The plot shows m\ (sohd hne) and /■„■ (dashed hne) as functions of the current 
quark mass rh in units of the physical value rhphys ■ 



tab. 0. To lowest order in (/jr — /4 }l j-n or, equivalently, in fn one has 



t _ f(o) _ 



4- 2 






(89) 



f(o) 



A larger value of m^ would therefore reduce the difference between f), and 

In fig. y we show the dependence of the pion mass and decay constant on 
the average current quark mass rh. These curves depend very little on the 
values of the initial parameters as demonstrated in tab. by /i . We observe 
a relatively large difference of 12 MeV between the pion decay constants at 
rh = TOphys and rh — Q. According to ( p9[ ) this difference is related to the mass 
of the sigma particle and will be modified in the three flavor case. We will later 
find that the critical temperature Tc for the second order phase transition in 
the chiral limit is almost independent of the initial conditions. The values of 
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/4 and Tc essentially determine the non-universal amplitudes in the critical 
scaling region (see below) . In summary, we find that the behavior of our model 
for small k is quite robust as far as uncertainties in the initial conditions at 
the scale fc$ are concerned. We will see that the difference of observables 
between non-vanishing and vanishing temperature is entirely determined by 
the flow of couplings in the range < A: ^ 2>T. 

10 Chiral phase transition of two flavor QCD 

Strong interactions in thermal equilibrium at high temperature T — as re- 
alized in early stages of the evolution of the Universe — differ in important 
aspects from the well tested vacuum or zero temperature properties. A phase 
transition at some critical temperature Tc or a relatively sharp crossover may 
separate the high and Iobl temperature physics .l3 Many experimental activi- 
ties at heavy ion coUideralj search for signs of such a transition. It was realized 
early that the transition should be closely related to a qualitative change in 
the chiral condensate according to the general observation that spontaneous 
symmetry breaking tends to be absmiLin a high temperature situation. A se- 
ries of stimulating contributionsOL3'Ej pointed out that for sufficiently small 
up and down quark masses, to„ and tti^, and for a sufficiently large mass of 
the strange quark, TOs, the chiral transition is expected to belong to the uni- 
versality class of the 0(4) Heisenberg model. This means that near the critical 
temperature only the pions and the sigma particle play a role for the hph avior 
of condensates and long distance correlation functions. It was suggestcdollj 
that a large correlation length may be responsible for important fluctu^ytipns 
or lead to a disoriented chiral coiidjensate.E3 This was even relatedllj' La to 
the spectacular "Centauro events"Ca observed in cosmic rays. The question 
how small to„ and nid would have to be in order to see a large correlation 
length near Tc and if this scenario could be realized for realistic values of the 
current quark masses remained, however, unanswered. The reason was the 
missing link between the universal behavior near Tc and zero current quark 
mass on one hand and the known physical properties at T == for realistic 
quark masses on the other hand. 

It is -the purpose of the remaining sections of these lectures to provide 
this link.L3 We present here the equation of state for two flavor QCD within 
an effective quark meson model. The equation of state expresses the chiral 
condensate (i/'''/') s-s ^ function of temperature and the average current quark 
mass m — {rnu-\'md)/'2.. This connects explicitly the universal critical behavior 
for T —f Tc and tti — > with the temperature dependence for a realistic value 
Wphys- Since our discussion covers the whole temperature range < T ^1.7 Tc 
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Figure 5. The plot shows the chiral condensate (t/iip) as a function of temperature T. Lines 
(a), (b), (c), (d) correspond at zero temperature to m„ = 0, 45 MeV, 135MeV, 230MeV, 
respectively. For each pair of curves the lower one represents the full T— dependence of 
(ipip) whereas the upper one shows for comparison the universal scaling form of the equation 
of state for the 0(4) Heisenberg model. The critical temperature for zero quark mass is 
Tc = 100.7 MeV. The chiral condensate is normalized at a scale fcj. ~ 620 MeV. 



we can fix rhphys such that the (zero temperature) pion mass is rriTr = 135 MeV. 
The condensate (V'V') plays here the role of an order parameter. Fig. || shows 
our results for {tpip) {T,m): Curve (a) gives the temperature dependence of 
(ipip') in the chiral limit m = 0. Here the lower curve is the full result for 
arbitrary T whereas the upper curve corresponds to the universal scaling 
form of the equation of state for the 0(4) Heisenberg model. We see perfect 
agreement of both curves for T sufficiently close to Tc = 100.7 MeV. This 
demonstrates the capability of our method to cover the critical behavior and, 
in particular, to reproduce the critical exponents of the 0(4)-model. We 
have determined (see below) the universal critical equation of state as well as 
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Table 2. The table shows the critical and "pseudocritical" temperatures for various values 
of the zero temperature pion mass. Here Tp^ is defined as the inflection point of {'ipipj (T) 

(2) 
whereas T™ is the location of the maximum of the sigma correlation length. 
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-110 


~ 130 


~ 150 


McV 


100.7 


113 


128 
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the non-universal amplitudes. This provides the full functional dependence 
of (V'V') {T,Th) for small T — Tc and m. The curves (6), (c) and (d) are 
for non-vanishing values of the average current quark mass rh. Curve (c) 
corresponds to mphys or, equivalently, niT^{T = 0) = 135 MeV. One observes 
a crossover in the range T = (1.2 — 1.5)rc. The 0(4) universal equation 
of state (upper curve) gives a reasonable approximation in this temperature 
range. The transition turns out to be much less dramatic than for rh = 0. 
We have also plotted in curve (6) the results for comparably small quark 
masses ~ 1 MeV, i.e. rh = Wphys/lO, for which the T = value of ttItt equals 
45 MeV. The crossover is considerably sharper but a substantial deviation 
from the chiral limit remains even for such small values of m. In order to 
facilitate comparison with lattice simulations which are typically performed 
for larger values of m^r we also present results for m-jriT = 0) = 230 MeV 
in curve (d). One may define a "pseudocritical temperature" Tpc associated 
to the smooth crossover as the inflection point of (V'V') {T) as often done in 

lattice simulations. Our results for this definition of Tpc are denoted by Tpc 
and are presented in tab. for the four different values of rh or, equivalently, 
TO7r(r = 0). The value for the pseudocritical temperature for ttt-tt = 230 MeV 
compares well with the lattice results for two flavor QCD (cf. the discussion 
below) . One should mention, though, that a determination of Tpc according 
to this definition is subject to sizeable numerical uncertainties for large pion 
masses as the curve in fig. H is almost linear around the inflection point for 
quite a large temperature range. A difficult point in lattice simulations with 
large quark masses is the extrapolation to realistic values of m^ or even to 
the chiral limit. Our results may serve here as an analytic guide. The overall 
picture shows the approximate validity of the 0(4) scaling behavior over a 
large temperature interval in the vicinity of and above Tc once the (non- 
universal) amplitudes are properly computed. 
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Figure 6. The plot shows m,r as a function of temperature T for three different values 
of the average light current quark mass m. The solid line corresponds to the realistic 
value m = rhphys whereas the dotted line represents the situation without explicit chiral 
symmetry breaking, i.e., m = 0. The intermediate, dashed line assumes m = rnpjjyg/lO. 



A second important result of our investigations is the temperature depen- 
dence of the space-Hke pion correlation length m~^{T). (We will often call 
TOtt (T) the temperature dependent pion mass since it coincides with the phys- 
ical pion mass for T — 0.) The plot for m^{T) in fig. ra again shows the second 
order phase transition in the chiral limit rfi — 0. For T < Tc the pions are 
massless Goldstone bosons whereas for T > Tc they form together with the 
sigma particle a degenerate vector of 0(4) whose mass increases as a function 
of temperature. For to = the behavior for small positive T — Tc is char- 
acterized by the critical exponent ly, i.e. mT^{T) = (^+) Tc{{T — Tc) /TcY 
and we obtain v = 0.787, ^+ = 0.270. For to > we find that m^(r) re- 
mains almost constant for T ^ Tc with only a very slight dip for T near 
Tc/2. For T > Tc the correlation length decreases rapidly and for T ^ Tc 
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the precise value of rh becomes irrelevant. We see that the universal critical 
behavior near Tc is quite smoothly connected to T = 0. The full functional 
dependence of ttItt (T, rh) allows us to compute the overall size of the pion cor- 
relation length near the critical temperature and we find ttt.^ (Tpc) — 1.7to^(0) 
for the realistic value rhphys- This correlation length is even smaller than the 
vacuum (T = 0) one and shows no indication for strong fluctuations of pions 
with long wavelength. It would be interesting to see if a decrease of the pion 
correlation length at and above Tc is experimentally observable. It should be 
emphasized, however, that a tricritical behavior with a massless excitation 
remains possible for three flavors. This would not be characterized by the 
universal behavior of the 0(4)-model. We also point out that the present 
investigation for the two flavor case does not take into account a speculaiijszE 
"effective restoration" of the axial Ua{^) symmetry at high temper ature.OLJ 
We will comment on these issues in the last section. In the next sections we 
will describe the formalism which leads to these results. 

11 Thermal equilibrium and dimensional reduction 

The extension of flow equations to thermal jequilibrium situations at non- 
vanishing temperature T is straightforward.Ea In the Euclidean formalism 
non-zero temperature results in (anti-)periodic boundary conditions fJCLc 
(fermionic) bosonic fields in the Euclidean time direction with periodicityE3 
1/T. This leads to the replacement 

in the trace in (O) when represented as a momentum integration, with a 
discrete spectrum for the zero component 

,.s _ j 21ttT for bosons ,„... 

^"^^^ \(2; + l)7rTforfermions. ^ ^ 

Hence, for T > a four-dimensional QFT can be interpreted as a three- 
dimensional model with each bosonic or fermionic degree of freedom now 
coming in an infinite number of copies labeled hy I € Z (Matsubara modes). 
Each mode acquires an additional temperature dependent effective mass term 
Qoil)- In a high temperature situation where all massive Matsubara modes 
decouple from the dynamics of the system one therefore expects to observe 
an effective three-dimensional theory with the bosonic zero modes as the 
only relevant degrees of freedom. In other words, if the characteristic length 
scale associated with the physical system is much larger than the inverse 
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temperature the compactified Euclidean "time" dimension cannot be resolved 
anymore. This phenomenon is known as "dimensional reduction" .Ell 

The formalism of the effective average action automatically provides the 
tools for a smooth decoupling of the massive Matsubara modes as the scale 
k is lowered from k ^ T to k -^ T. It therefore allows us to directly link 
the low-T, four-dimensional QFT to the effective three-dimensional high-T 
theory. The replacement ( po| ) in (O) manifests itself in the flow equations ^74), 
(0) — (^9|) only through a change to T-dependent threshold functions. For 
instance, the dimensionless functions /^(w; ?/$) defined in ( p3[ ) are replaced by 



where q^ = Qq+q'^ and go = 2t:IT . A list of the various temperatur£_dependent 
threshold functions appearing in the flow equations can be foundH There we 
also discuss some subtleties regarding the definition of the Yukawa coupling 
and the anomalous dimensions for T ^ Q. In the limit k ^ T the sum over 
Matsubara modes approaches the integration over a continuous range of qo 
and we recover the zero temperature threshold function lf^{w,ri^). In the 
opposite limit k <^ T the massive Matsubara modes {I ^ 0) are suppressed 
and we expect to find a d— 1 dimensional behavior of ZjJ. In fact, one obtains 
from (H) 

Z^(w, T/k; ri^) ~ li{w] 77$) for T < A: , 

Tv, 1 (93) 

li{w, T/k; 7j^) ~ -^^l^-\w; 77$) for T > A: . 
k Vd 

For our choice of the infrared cutoff function Rk, (||), the temperature de- 
pendent Matsubara modes in l'l{w, T/k; ry$) are exponentially suppressed for 
T <^ k whereas the behavior is more complicated for other threshold func- 
tions appearing in the flow equations ([74[), ( |77| ) — ([79[). Nevertheless, all 
bosonic threshold functions are proportional to T/k for T ^ k whereas those 
with fermionic contributions vanish in this limitP] This behavior is demon- 
strated in fig. 1^ where we have plotted the quotients lf{w,T/k)/lf{w) and 
l{ {w,T/k)/l\ (w) of bosonic and fermionic threshold functions, respec- 
tively. 



''For the present choice of ij^ the temperature dependence pfithe threshold functions is 
considerably smoother than in that of previous investigations 123- 
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Figure 7. The plot shows the temperature dependence of the bosonic (a) and the fcrmionic 
(b) threshold functions Z|(ui, T/k) and r-^ ' (w, T/k), respectively, for different values of the 
dimensionless mass term w. The solid line corresponds to lo = whereas the dotted ones 
correspond to ui = 0.1, w = 1 and ui = 10 with decreasing size of the dots. For T ^ k the 
bosonic threshold function becomes proportional to T/k whereas the fermionic one tends 
to zero. In this range the theory with properly rescaled variables behaves as a classical 
three-dimensional theory. 



One observes that for k ^ T both threshold functions essentiaUy behave as 



faro: submitted to World Scientific on December 18, 1998 



50 



for zero temperature. For growing T or decreasing k this changes as more 
and more Matsubara modes decouple until finally all massive modes are sup- 
pressed. The bosonic threshold function if shows for fc <C T the linear de- 
pendence on T/k derived in (p3[). In particular, for the bosonic excitations 
the threshold function for w <^ 1 can be approximated with reasonable ac- 
curacy by lt{w;r]^) for T/k < 0.25 and by {AT / k)ll{w; r]^) for T/k > 0.25. 
The fermionic threshold function ll tends to zero for k <^ T since there 
is no massless fermionic zero mode, i.e. in this limit all fermionic contribu- 
tions to the flow equations are suppressed. On the other hand, the fermions 
remain quantitatively relevant up to T/k ~ 0.6 because of the relatively long 
tail in fig. 0b. The transition from four to three-dimensional threshold func- 
tions leads to a smooth dimensional reduction as k is lowered from k ^ T to 
k -^ Tl Whereas for k :^ T the model is most efficiently described in terms 
of standard four-dimensional fields $ a choice of rescaled three-dimensional 
variables $3 = <i>/\/T becomes better adapted for k <$^ T. Accordingly, for 
high temperatures one will use a potential 

k k 

U3it,P3) = 7^u{t,~p) ; p3 = -p . (94) 

In this regime Tk~,o corresponds to the free energy of classical statistics and 
rfc>o is a classical coarse grained free energy. 

For our numerical calculations at non-vanishing temperature we exploit 
the discussed behavior of the threshold functions by using the zero tempera- 
ture flow equations in the range k > lOT. For smaller values of k we approx- 
imate the infinite Matsubara sums (cf. (^)) by a finite series such that the 
numerical uncertainty at fc = lOT is better than 10^^. This approximation 
becomes exact in the limit k <C lOT. 



12 The quark meson model at T 7^ 

So far we have considered the relevant fluctuations that contribute to the 
flow of Ffc in dependence on the scale k. In a thermal equilibrium situation 
Ffc also depends on the temperature T and one may ask for the relevance 
of thermal fluctuations at a given scale fc. In particular, for not too high 
values of T the "initial condition" Ffc^ for the solution of the flow equations 
should essentially be independent of temperature. This will allow us to flx 
Ffc^ from phenomenological input at T = and to compute the tempera- 
ture dependent quantities in the infrared (fc ^ 0). We note that the thermal 
fluctuations which contribute to the r.h.s. of the flow equation for the meson 
potential ( |74| ) are effectively suppressed for T ^ fc/4 as discussed in detail in 
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the last section. Clearly for T ^ fc$/3 temperature effects become important 
at the compositeness scale. We expect the linear quark meson model with 
a compositeness scale fc$ ~ 600 MeV to be a valid description for two flavor 
QCD below a temperature of about 170 MeV. We note that there will be 
an effective temperature dependence of F^^ induced by the fluctuations of 
other degrees of freedom besides the quarks, the pions and the sigma which 
are taken into account here. We will comment on this issue in the last sec- 
tion. For realistic three flavor QCD the thermal kaon fluctuations will become 
important for T > 170 MeV. 

We compute the quantities of interest for temperatures T < IjZDMeV 
by numerically solving the T-dependent version of the flow equationsE3 (174|), 
( |77| ) — ( [79[ ) by lowering k from fc$ to zero. For this range of temperatures we 
use the initial values as given in the first line of tab. |l|. This corresponds 
to choosing the zero temperature pion mass and the pion decay constant 
{fir — 92.4 MeV for m-rr = 135 MeV) as phenomenological input. The only 
further input is the constituent quark mass Mg which we vary in the range 
Mq ~ 300 — 350 MeV. We observe only a minor dependence of our results on 
Mg for the considered range of values. In particular, the value for the critical 
temperature T^ of the model remains almost unaffected by this variation. 

We have plotted in fig. ^ the renormalized expectation value 2aQ of the 
scalar field as a function of temperature for three different values of the average 
light current quark mass m. (We remind that 2ao{T = 0) = f^.) For to = 
the order parameter ctq of chiral symmetry breaking continuously goes to zero 
for T ^ Tc = 100.7 MeV characterizing the phase transition to be of second 
order. The universal behavior of the model for small T — Tc and small to 
is discussed in more detail in the following section. We point out that the 
value of Tc corresponds to /4 = 80.8 MeV, i.e. the value of the pion decay 
constant for to = 0, which is significantly lower than /^r = 92.4 MeV obtained 
for the realistic value mphys- If we would fix the value of the pion decay 
constant to be 92.4 MeV also in the chiral limit (to = 0), the value for the 
critical temperature would raise to 115 MeV. The nature of the transition 
changes qualitatively for to 7^ where the second order transition is replaced 
by a smooth crossover. The crossover for a realistic TOphys or mTr{T = 0) = 
135 MeV takes place in a temperature range T ~ (120 - 150) MeV. The 
middle curve in fig. ra corresponds to a value of to which is only a tenth of 
the physical value, leading to a zero temperature pion mass TOtt = 45 MeV. 
Here the crossover becomes considerably sharper but there remain substantial 
deviations from the chiral limit even for such small quark masses m ~ 1 MeV. 
The temperature dependence of tti^t has already been mentioned (see fig. ra) for 
the same three values of to. As expected, the pions behave like true Goldstone 
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Figure 8. The expectation value 2(To is shown as a function of temperature T for three 
different values of the zero temperature pion mass. 



bosons for to = 0, i.e. their mass vanishes for T < Tc- Interestingly, ?7i^ 
remains almost constant as a function of T for T < Tc before it starts to 
increase monotonically. We therefore find for two flavors no indication for a 
substantial decrease of to..^- around the critical temperature. 

The dependence of the mass of the sigma resonance m^ on the tempera- 
ture is displayed in fig. for the above three values of m 
In the absence of explicit chiral symmetry breaking, m = 0, the sigma mass 
vanishes for T < Tc. For T < Tc this is a consequence of the presence of 
massless Goldstone bosons in the Higgs phase which drive the renormalized 
quartic coupling A to zero. In fact, A runs linearly with k for T ^ fc/4 and 
only evolves logarithmically for T ^ fc/4. Once to 7^ the pions acquire a 
mass even in the spontaneously broken phase and the evolution of A with fc 
is effectively stopped at fc ~ ttItt- Because of the temperature dependence 
of (Jo^k=o (cf. fig. gf) this leads to a monotonically decreasing behavior of ttIo- 
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Figure 9. The plot shows the m^ as a function of temperature T for three different values 
of the zero temperature pion mass. 



with T for T ^ Tc- This changes into the expected monotonic growth once the 
system reaches the symmetric phase for T > T^. For low enough rh one may 
use the minimum of mo-(T) for an alternative definition of the (pseudo-)critical 
temperature denoted as Tpc . Tab. || in the introduction shows our results 
for the pseudocritical temperature for different values of m or, equivalently, 
r7z^(T = 0). For a zero temperature pion mass ttItt = 135 MeV we find Tpc = 
128 MeV. At larger pion masses of about 230 MeV we observe no longer a 
characteristic minimum for m„ apart from a very broad, slight dip at T ~ 
90 MeV. A comparison of our results with lattice data is given below in the 
next section. 

Our results for the chiral condensate (?/"/') as a function of temperature 
for different values of the average current quark mass are presented in fig. ^. 
We will compare (^V') (Tj'rh) with its universal scaling form for the 0(4) 
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Figure 10. The plot shows the derivative of the meson potential U{T) with respect to the 
renormalizcd field (pji = (Z$p/2)^" for different values of T. The first curve on the left 
corresponds to T = 175 MeV. The successive curves to the right differ in temperature by 
AT = 10 MeV down to T = 5 MeV. 



Heisenberg model in the following section. 

Our ability to compute the complete temperature dependent effective me- 
son potential U is demonstrated in fig. M where we display the derivative of 
the potential with respect to the renormalized field (pR = (Z$p/2)^/^, for dif- 
ferent values of T. The curves cover a temperature range T — {5— 175) MeV. 
The first one to the left corresponds to T = 175 MeV and neighboring curves 
differ in temperature by AT — 10 MeV. One observes only a weak dependence 
of dU{T)/d(j)R on the temperature for T SS 60 MeV. Evaluated for (j)ji = ctq 
this function connects the renormalized field expectation value with TO7r(T), 
the source j and the mesonic wave function renormalization Z$ (T) according 
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to 

^^i^R - ^o) = ^1^ = 4ao(T)m^(r) . (95) 

d(pR ZH^T) 

We close this section with a short assessment of the vahdity of our ef- 
fective quark meson model as an effective description of two flavor QCD at 
non-vanishing temperature. The identification of qualitatively different scale 
intervals which appear in the context of chiral symmetry breaking, as pre- 
sented in the preceding sections for the zero temperature case, can be gener- 
alized to T 7^ 0: For scales below kg, there exists a hybrid description in terms 
of quarks and mesons. For k^sB < k ^ 600 MeV chiral symmetry remains un- 
broken where the symmetry breaking scale k^sB {T) decreases with increasing 
temperature. Also the constituent quark mass decreases with T. The run- 
ning Yukawa coupling depends only mildly on temperature for T ^ 120 MeV. 
(Only near the critical temperature and for m — the running is extended 
because of massless pion fluctuations.) On the other hand, for k ^ AT the 
effective_^three-dimensional gauge coupling increases faster than at T = 
leading to an increase of Aqcd(T) with T. As k gets closer to the scale 
^QCT>{T) it is no longer justified to neglect in the quark sector confinement 
effects which go beyond the dynamics of our present quark meson model. 
Here it is important to note that the quarks remain quantitatively relevant 
for the evolution of the meson degrees of freedom only for scales k ^ r/0.6 
(cf. fig. 0). In the limit k <C T/0.6 all fermionic Matsubara modes decou- 
ple from the evolution of the meson potential according to the temperature 
dependent version of ([74|). Possible sizeable confinement corrections to the 
meson physics may occur if Aqcd {T) becomes larger than the maximum of 
Mq{T) and T/0.6. This is particularly dangerous for small to in a temperature 
interval around Tc- Nevertheless, the situation is not dramatically different 
from the zero temperature case since only a relatively small range of k is 
concerned. We do not expect that the neglected QCD non-localities lead to 
qualitative changes. Quantitative modifications, especially for small to and 
\T — Tc| remain possible. This would only effect the non-universal amplitudes 
which will be discussed in the next section. The size of these corrections de- 
pends on the strength of (non-local) deviations of the quark propagator and 
the Yukawa coupling from the values computed in the quark meson model. 

13 Critical behavior near the chiral phase transition 

In this section we study the linear quark meson model in the vicinity of the 
critical temperature T^ close to the chiral limit to = 0. In this region we 
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find that the sigma mass m~^ is much larger than the inverse temperature 
T~^, and one observes an effectively three-dimensional behavior of the high 
temperature quantum field theory. We also note that the fermions are no 
longer present in the dimensionally reduced system as has been discussed 
above. We therefore have to deal with a purely bosonic 0(4)-symmetric linear 
sigma model. At the phase transition the correlation length becomes infinite 
and the effective three-dimensional theory is dominated by classical statistical 
fluctuations. In particular, the critical exponents which describe the singular 
behavior of various quantities near the second order phase transition are those 
of the corresponding classical system. 

Many properties of this system are universal, i.e. they only depend on 
its symmetry (0(4)), the dimensionality of space (three) and its degrees of 
freedom (four real scalar components). Universality means that the long- 
range properties of the system do not depend on the details of the specific 
model like its short distance interactions. Nevertheless, important properties 
as the value of the critical temperature are non-universal. We emphasize that 
although we have to deal with an effectively three-dimensional bosonic theory, 
the non-universal properties of the system crucially depend on the details of 
the four-dimensional theory and, in particular, on the fermions. 

Our aim is a computation of the critical equation of state which relates 
for arbitrary T near T^ the derivative of the free energy or effective potential 
U to the average current quark mass m. The equation of state then permits to 
study the temperature and quark mass dependence of properties of the chiral 
phase transition. 

At the critical temperature and in the chiral limit there is no scale present 
in the theory. In the vicinity of Tc and for small enough rh one therefore ex- 
pects a scaling behavior of the effective average potential Uk and accordingly 
a universal scaling form of the equation of state.E3 There are only two in- 
dependent scales close to the transition point which can be related to the 
deviation from the critical temperature, T — T^,, and to the explicit symmetry 
breaking through the quark mass m. As a consequence, the properly rescaled 
potential can only depend on one scaling variable. A possible choice for the 
parameterization of the_Efiscaled "unrenormalized" potential is the use of the 



Widom scaling variable 



(T - Te) /T, 
{2Wo/Tc 



\l//3 



(96) 



Here /3 is the critical exponent of the order parameter ctq in the chiral limit 



m 



(see (llOq)). With U'{p = 2a^) = j/(2cto) the Widom scaling form of 
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the equation of state readsE3 

s 



where the exponent S is related to the behavior of the order parameter 



according to (102). The equation of state ( pT\) is written for convenience 
directly in terms of four-dimensional quantities. They are related to the 
corresponding effective variables of the three-dimensional theory by appro- 
priate powers of Tc- The source j is determined by the average current 
quark mass m as j = 2to^ rh. The mass term at the compositeness scale, 
mf. , also relates the chiral condensate to the order parameter according to 
i^ipip) = — 2to|^ (ctq ~ ™)- The critical temperature of the linear quark meson 
model was found above to be Tc = 100.7 MeV. 

The scaling function / is universal up to the model specific normalization 
of X and itself. Accordingly, all models in the same universality class can be 
related by a rescaling of ctq and T — Tc- The non-universal normalizations for 
the quark meson model discussed here are defined according to 

f{Q) = D , /(-i?-!/'^) = . (98) 

We find D = 1.82 • 10"'*, B = 7.41 and our resuh for (3 is given in tab. |. 
Apart from the immediate vicinity of the_zj3ro of f{x) we find the following 
two parameter fit for the scaling function,E3 

ffitix) = 1.816 • 10"*(1 + 136.1 x)^ (1 + 160.9 6'x)^ 

(99) 

(1 + 160.9 (0.9446 0-i/(^-2-A) ^)7-2-A 

to reproduce the numerical results for / and df/dx at the 1 — 2% level with 
e = 0.625 (0.656), A = -0.490 (-0.550) for x > (x < 0) and 7 as given in 
tab. 0. The universal properties of the scaling function can be compared with 
results obtained by other methods for the three-dimensional 0(4) Heisenberg 



model. In fig. 11 we displai^ our results along with those obtained from lattice 
Monte Carlo simulation,E3 second order epsilon expansions and mean field 
theory. We observe a good agreement of average action, lattice and epsilon 
expansion results within a few per cent for T < Tc- Above Tc the average 
action and the lattice curve go quite close to each other with a substantial 
deviation from the epsilon expansion and mean field scaling function. (We 
note that the question of a better agreement of the curves for T < Tc or 
T > Tc depends on the chosen non-universal normalization conditions for x 
and / (cf. (H)).) 
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Figure 11. The figure shows a comparison of our results, denoted by "average action", with 
results of other methods for the scaling function of the three-dimensional 0{4) Heisenberg 
model. We have labeled the axes for convenience in terms of the expectation value ctq and 
the source j of the corresponding four— dimensional theory. The constants B and D specify 
the non— universal amplitudes of the model (cf. (bsL The curve labeled br-|"MC" represents 
a fit to lattice Monte Carlo data. The second order epsilon expansionEj and mean field 
results are denoted by "e" and "mf", respectively. Apart from our results the curves are 
taken from Rcf.LJ. 



Before we use the scaling lunction f{x) to discuss the general temperature 
and quark mass dependent case, we consider the hniits T = Tc and to = 0, re- 
spectively. In these limits the behavior of the various quantities is determined 
solely by critical amplitudes and exponents. In the spontaneously broken 
phase (T < Tc) and in the chiral limit we observe that the renormalized and 
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unrenormalized order parameters scale according to 



2ao(r) _ ^2^^V2 (T,^ 






T, V T, 



vjl 

(100) 



respectively, with E — 0.814 and the value of i? given above. In the symmetric 
phase the renormalized mass m = m^ = Too- and the unrenormalized mass 
TO = ZJ m behave as 



m{T) _ /.+^-l fT-Tc 



Tr ^" ^ V T, 



■m[T) _ .^+^^-l/2 (T-Tc 



n ' ' V T, 



./2 (101) 



where 5"*" = 0.270, C+ = 2.79. For T ^ Tc and non-vanishing current quark 
mass we have 

with the value of D given above. 

Though the five amplitudes E, B, ^'^ , C^ and D are not universal there 
are ratios of amplitudes which are invariant under a rescaling of ctq and T — Tc. 
Our results for the universal amplitude ratios are 

i?x = C+DB^-^ = 1.02 , 

R^ = {^+f/''D^/^^+^'>B = 0.852 , (103) 

^+E ^ 0.220 . 

Those for the critical exponents are given in tab. 0. Here the value of rj 
is obtained fcom the temperature dependent version of ([7^^^ the critical 
temperature^ For comparison tab. ^ also gives the resultsBEll of a 3d per- 
turbative computation as well as lattice Monte Carlo resultsEa which have been 
used for the lattice form of the scaling function in fig. ^. [] There are only 
two independent amplitudes and critical exponents, respectively. They are 



'See also Ref.E3 and references therein for a recent calculation of critical exponents using 
similar n»ej;]|*e|ds as in this work. For high precision estimates of the critical exponents see 
also Refsl 
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Table 3. The table shows the critical exponents corresponding to the three-dimensional 
0(4)— Heisenberg model. Our results are denoted by "AA" whereas "3d PT" lieJjels the 
exponents obtained from a 3d perturbative exreapsion using Pade techniques ££r|£j The 
bottom line contains lattice Monte Carlo results.Ej 





v 


7 


S 


(3 


V 


AA 


0.787 


1.548 


4.80 


0.407 


0.0344 


3dPT 


0.73(2) 


1.44(4) 


4.82(5) 


0.38(1) 


0.03(1) 


MC 


0.7479(90) 


1.477(18) 


4.851(22) 


0.3836(46) 


0.0254(38) 



relatecLby the usual scaling relations of the three-dimensional scalar 0{N)- 
modelEO which we have explicitly verified by the independent calculation of 
our exponents. 

We turn to the discussion of the scaling behavior of the chiral condensate 
(f/'f/A for the general case of a temperature and quark mass dependence. In 
fig. ra we have displayed our results for the scaling equation of state in terms 
of the chiral condensate 



(^^) = 



-mlTr, 



JJTI 



+ J 



(104) 



as a function of T/Tc = 1 + x{j/T'^f{x)Y^^^ for different quark masses or, 
equivalently, different values of ]. The curves shown in fig. a correspond 
to quark masses rh — 0, m — TOphys/10, m — mphys and m — 3.5mphys 
or, equivalently, to zero temperature pion masses ttItt = 0, ttItt = 45 MeV, 
m.^ = 135 MeV and TOtt — 230 MeV, respectively (cf. fig. ^). One observes 
that the second order phase transition with a vanishing order parameter at 
Tc for TO == is turned into a smooth crossover in the presence of non-zero 
quark masses. 



The scaling form (104) for the chiral condensate is exact only in the limit 
T ^ Tc, J — > 0. It is interesting to find the range of temperatures and quark 
masses for which i^pip) approximately shows the scaling behavior (104). This 
can be inferred from a comparison (see fig. p[) with our full non-universal 
solution for the T and j dependence of {i^ip)- For TOtt — one observes 
approximate scaling behavior for temperatures T ^ 90 MeV. This situation 
persists up to a pion mass of to-tt — 45 MeV. Even for the realistic case, 
TO-TT = 135 MeV, and to a somewhat lesser extent for to-jt — 230 MeV the scaling 
curve reasonably reflects the physical behavior for T ^ Tc- For temperatures 
below Tc, however, the zero temperature mass scales become important and 
the scaling arguments leading to universality break down. 



faro: submitted to World Scientific on December 18, 1998 



61 



The above comparison may help to shed some hght on the use of univer- 
sahty arguments away from the critical temperature and the chiral limit. One 
observes that for temperatures above Tc the scaling assumption leads to quan- 
titatively reasonable results even for a pion mass almost twice as large as the 
physical value. This in turn has been used for two flavor lattice QCD as the- 
oretical input to guide extrapolation of results to light current quark masses. 
From simulations based on a range of pion masses 0.3 ^ m^r/mp ^ 0.7 and 
temperatures Q < T <, 250 MeV a "pseudocritical temperature" f£ approxi- 
mately 140 MeV with a weak quark mass dependence is reported.Ell Here the 
"pseudocritical temperature" Tpc is defined as the inflection point of (i/'V') as-a 
function of temperature. The values of the lattice action parameters used irO 
with iVi = 6 were am = 0.0125, 6/5^ = 5.415 and am = 0.025, 6/3^ = 5.445. 
For comparison with lattice data we have displayed in fig. g the temperature 
dependence of the chiral condensate for a pion mass m,.„: = 230 MeV. From the 
free energy of the linear quark meson model we obtain in this case a pseudocrit- 
ical tepiperature of about 150 MeV in reasonable agreement with the results 
of Ref.Ej. In contrast, for the critical temperature in the chiral limit we obtain 
Tc = 100.7 MeV. This value is considerably smaller than the lattice results 
of about (140 — 150) MeV obtained by extrapolating to zero quark mass in 
Ref.Lj. We point out that for pion masses as large as 230 MeV the condensate 
(ipipj {T) is almost linear around the inflection point for quite a large range 
of temperature. This makes a precise determination of Tc somewhat difficult. 
Furthermore, fig. |5| shows that the scaling form of (V'V') (^) underestimates 
the slope of the physical curve. Used as a fit with Tc as a parameter this can 
lead to an overestimate of the pseudocritical temperature in the chiral limit. 
We also mention here the results of RefB3. There two values of the pseudo- 
critical temperature, Tpc ~ 150(9) MeV and Tpc = 140(8), corresponding to 
am = 0.0125, Q/g^ = 5.54(2) and am = 0.00625, Q/g'^ = 5.49(2), respectively, 
(both for Nt — 8) were computed. These values show a somewhat stronger 
quark mass dependence of Tpc and were used for a linear extrapolation to the 
chiral limit yielding Tc = 128(9) MeV. 

The linear quark meson model exhibits a second order phase transition for 
two quark flavors in the chiral limit. As a consequence the model predicts a 
scaling behavior near the critical temperature and the chiral limit which can, 
in principle, be tested in lattice simulations. For the quark masses used in the 
present lattice studies the order and universality class of the transition in two 
flavor QCD remain a partially open question—. Though there are results from 
the lattice giving support for critical scalingHj'O there are also recent simu- 
lations with two fLaYQis that reveal significant finite size effects and problems 
with 0(4) scaling.0Ea 
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14 Additional degrees of freedom 

So far we have investigated the chiral phase transition of QCD as described 
by the hnear 0(4)-niodel containing the three pions and the sigma resonance 
as well as the up and down quarks as degrees of freedom. Of course, it is clear 
that the spectrum of QCD is much richer than the states incorporated in our 
model. It is therefore important to ask to what extent the neglected degrees of 
freedom like the strange quark, strange (pseudo)scalar mesons, (axial)vector 
mesons, baryons, etc., might be important for the chiral dynamics of QCD. 
Before doing so it is perhaps instructive to first look into the opposite di- 
rection and investigate the difference between the linear quark meson model 
described here and chiral perturbation theory based on the non-linear sigma 
model.Ej In some sense, chiral perturbation theory is the minimal model of 
chiral symmetry breaking containing only the Goldstone degrees of freedom. 
By construction it is therefore only valid in the spontaneously broken phase 
and can not be expected to yield realistic results for temperatures close to Tc 
or for the symmetric phase. However, for small temperatures (and momentum 
scales) the non-linear model is expected to describe the low-energy and low- 
temperature limit of QCD reliably as an expansion in powers of the light-aaapk 
masses. For vanishing temperature it has been demonstrated recentlyE^'ESEJ 
that the results of chiral perturbation theory can be reproduced within the 
linear meson model once certain higher dimensional operators in its effective 
action are taken into account for the three flavor case. Moreover, some of the 
parameters of chiral perturbation theory (^4, . . . , Lg) can be expressed and 
therefore also numerically computed in terms of those of the linear model. For 
non-vanishing temperature one expects agreement only for low T whereas de- 
viations from chiral perturbation theory should become large close to Tc- Yet, 
even for T <^Tc small quantitative deviations should exist because of the con- 
tributions of (constituent) quark and sigma meson fluctuations in the linear 
model which are not taken into account in chiral perturbation theory. 

FromEZI we infer the three-loop result for the temperature dependence of 
the chiral condensate in the chiral limit for N light flavors 

— — f A^^ — 1 T^ 

(V'V) (T)xPT = (V'V')^P^ (0) 1 - 




+ iV(iV2 - 1) 

y iiX-Q / 1 1 I 

(105) 
The scale Fi can be determined from the _D-wave isospin zero tttt scattering 
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Figure 12. The plot displays the chiral condensate (ipip) as a function of T/ f)^ . The solid 
line corresponds to our results for vanishing average current quark mass m = whereas 
the dashed line shows the corresponding three-loop chiral perturbation theory result for 
Ti =470MeV. 



length and is given by Ti = (470 ± 100) MeV. The constant Fq is (in the 
chiral limit) identical to the pion decay constant Fq — frr — 80.8 MeV. In 
fig. |l^ we have plotted the chiral conde nsate as a function of T/Fq for both, 
chiral perturbation theory according to (105) and for the linear quark meson 
model. As expected the agreement for small T is very good. Nevertheless, 
the anticipated small numerical deviations present even for T <^ Tc due to 
quark and sigma meson loop contributions are manifest. For larger values 
of T, say for T ^ 0.8/4 the deviations become significant because of the 
intrinsic inability of chiral perturbation theory to correctly reproduce the 
critical behavior of the system near its second order phase transition. 

Within the language of chiral perturbation theory the neglected effects 
of thermal quark fluctuations may be described by an effective temperature 
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dependence of the parameter Fo{T). We notice that the temperature at which 
these corrections become important equals approximately one third of the 
constituent quark mass Mq{T) or the sigma mass m„{T)^ respectively, in 
perfect agreement with fig. R. As suggested by this figure the onset of the 
effects from thermal fluctuations of heavy particles with a T-dependent mass 
mH{T) is rather sudden for T ^ mH{T)/i. These considerations also apply to 
our two flavor quark meson model. Within full QCD we expect temperature 
dependent initial values at A:$. 

The dominant contribution to the temperature dependence of the initial 
values presumably arises from the influence of the mesons containing strange 
quarks as well as the strange quark itself. Here the quantity m^ seems to be 
the most important one. (The temperature dependence of higher couplings 
like A(T) is not very relevant if the IR attractive behavior remains valid, i.e. if 
Z-i>,k^ remains small for the range of temperatures considered. We neglect a 
possible T-dependence of the current quark mass rh.) In particular, for three 
flavors the potential Uk^ contains a term 

- -Vk^ (det $ + det $^) = -Vk^ Vs^uu^dd + ■■■ (106) 

which reflects the axial Ua{^) anomaly. It yields a contribution to the effective 
mass term proportional to the expectation value {^ss) = <^s, i-S. 

^ml^^-'^Vk^^s. (107) 

Both, I7fc^ and (p^, depend on T. We expect these corrections to become 
relevant only for temperatures exceeding m/<-(r)/3 or Ms{T)/i. We note 
that the temperature dependent kaon and strange quark masses, mK{T) and 
Ms{T), respectively, may be somewhat different from their zero temperature 
values but we do not expect them to be much smaller. A typical value for 
these scales is around 500 MeV. Correspondingly, the thermal fluctuations 
neglected in our model should become important for T > 170 MeV. It is 
even conceivable that a discontinuity appears in ips (T) for sufficiently high T 
(say T ~ 170 MeV). This would be reflected by a discontinuity in the initial 
values of the 0(4)-model leading to a first order transition within this model. 
Obviously, these questions should be addressed in the framework of the three 
flavor SUl{3) x SUii{3) quark meson model. Work in this direction is in 
progress. 

We note that the temperature dependence of T'{T)ifs{T) is closely re- 
lated to the question-oL^n effective high temperature restoration of the ax- 
ial Ua{^) symmetryO'Lj The rj' mass term is directly proportional to this 
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combination J^ m'^,(T) — m'^{T) ~ ^V{T)(psiT). Approximate /7yi(l) restora- 
tion would occur if (ps (T) or T'{T) would decrease sizeable for large T. For re- 
alistic QCD this question should be addressed by a three flavor study. Within 
two flavor QCD the combination T'k(ps is replaced by an effective anomalous 
mass term i7^ . The temperature dependence of T7'^'(T) could be studied by 
introducing quarks and the axial anomaly in the two flavor matrix model of 
Ref.LD. We add that this questiop-has also been studied within full two flavor 
QCD in lattice simulations Ot^Ej So far there does not seem to be much evi- 
dence for a restoration of the f/A(l) symmetry near Tc but no final conclusion 
can be drawn yet. 

To summarize, we have found that the effective two flavor quark meson 
model presumably gives a good description of the temperature effects in two 
flavor QCD for a temperature range T <, 170 MeV. Its reliability should be 
best for low temperature where our results agree with chiral perturbation the- 
ory. However, the range of validity is considerably extended as compared to 
chiral perturbation theory and includes, in particular, the critical temperature 
of the second order phase transition in the chiral limit. We have explicitly 
connected the universal critical behavior for small jT — Tc| and small cur- 
rent quark masses with the renormalized couplings at T = and realistic 
quark masses. The main quantitative uncertainties from neglected fluctua- 
tions presumably concern the values of /4 and Tc which, in turn, influence 
the non-universal amplitudes B and D in the critical region. We believe that 
our overall picture is rather solid. Where applicable our results compare well 
with numerical simulations of full two flavor QCD. 

15 Conclusions 

Our conclusions may be summarized in the following points: 

1. The connection between short distance perturbative QCD and long dis- 
tance meson physics by analytical methods seems to emerge step by step. 
The essential ingredients are nonperturbative flow equations as approxi- 
mations of exact renormalization group equations and a formalism which 
allows a change of variables by introducing meson-like composite fields. 

2. The relevance of meson-like ^pip composite objects is established in this 
framework. A typical scale where the mesonic bound states appear is 
the compositeness scale k^ ~ (600 — 700) MeV. The occurrence of chiral 
symmetry breaking depends on the effective meson mass mfc^ at the com- 
positeness scale. For a certain range of values of the ratio to| /A:| spon- 
taneous chiral symmetry breaking is induced by quark fluctuations. A 
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definite analytical establishment of spontaneous chiral symmetry break- 
ing from "first principles" (i.e., short distance QCD) still awaits a reliable 
calculation of this ratio. 

Phcnomcnologically, the ratio m^ /k\ may be determined from the value 
of the constituent quark mass in units of the pion decay constant /jr. For 
this value the ratios /^/fctj, and (i/'V') /^l can be computed. Together 
with an earlier estimate of fc$ this yields rather encouraging results: /^ ~ 
100 MeV and (V'V') ^ -(190 MeV)^ for two flavor QCD. It wiU be very 
interesting to generalize these results to the realistic three flavor case. 

The formalism presented here naturally leads to an effective linear quark 
meson model for the description of mesons below the compositeness scale 
k^. For this model the standard non-linear sigma model of chiral per- 
turbation theory emerges as a low energy approximation due to the large 
sigma mass. 

The old puzzle about the precise connection between the current and the 
constituent quark mass reveals new interesting aspects in this formalism. 
In the context of the effective average action one may define these masses 
by the quark propagator at zero or very small momentum, either in the 
quark gluon picture (current quark mass) or the effective quark meson 
model (constituent quark mass). There is no conceptual difference be- 
tween the two situations. As a function of k the running quark mass 
smoothly interpolates between the standard current quark mass mq for 
high k and the standard constituent quark mass Mq for low k. (This 
holds at least as long as the minimum of the effective scalar potential is 
unique.) A rapid quantitative change occurs for k ~ (400 — 500) MeV 
because of the onset of chiral symmetry breaking. For fc = one expects 
this behavior to carry over to the momentum dependence of the quark 
propagator: For small q^ the inverse quark propagator is dominated by 
Mq. In contrast, for high q^ the constant term in the inverse propagator 
is reduced to ruq since q^ replaces k^ as an effective infrared cutoff. One 
expects a smooth interpolation between the two limits and it would be 
interesting to know the form of the propagator for momenta in the tran- 
sition region. Furthermore, the symmetry breaking source term which 
determines the pion mass in the linear or non-linear meson model can 
be related to the current quark mass at the compositeness scale. This 
constitutes a bridge between the low energy meson properties and the 
running quark mass at a scale which is not too far from the validity of 
perturbation theory. 



faro: submitted to World Scientific on December 18, 1998 67 



A particular version of the Nambu-Jona-Lasinio model appears in our 
formalism as a limiting case (infinitely strong rcnormalized Yukawa cou- 
pling at the compositencss scale kg,). Here the ultraviolet cutoff which is 
implicit in the NJL model is dictated by the momentum dependence of 
the infrared cutoff function Rk in the effective average action. The char- 
acteristic cutoff scale is fc$. In this context our results can be interpreted 
as an approximative solution of the NJL model. Our method includes 
many contributions beyond the leading order contribution of the 1/Nc 
expansion. 

Based on a satisfactory understanding of the meson properties in the 
vacuum we have described their behavior in a thermal equilibrium situa- 
tion. Our method should remain valid for temperature below ~ 170 MeV. 
This extends well beyond the validity of chiral perturbation theory. In 
particular, in the chiral limit of vanishing quark masses we can describe 
the universal critical behavior near a second order phase transition of 
two flavor QCD. The universal behavior is quantitatively connected to 
observed quantities at zero temperature and realistic quark masses. 
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